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Abstract 

We establish existence and uniqueness for the martingale problem associated with 
a system of degenerate SDE's representing a catalytic branching network. For example, in 
the hypercyclic case: 



dxi'^ = h{Xt)dt + \j2^,{Xt)X['+^^x'thBl Xf^ > 0, z = 1, . . . , ci, 

where X'^'^^^^ = X'^'^\ existence and uniqueness is proved when 7 and h are continuous 
on the positive orthant, 7 is strictly positive, and 6^ > on {xi = 0}. The special case 
d = 2, bi = 9i — Xi is required in work of [DGHSS] on mean fields limits of block averages 
for 2-type branching models on a hierarchical group. The proofs make use of some new 
methods, including Cotlar's lemma to establish asymptotic orthogonality of the derivatives 
of an associated semigroup at different times, and a refined integration by parts technique 
from [DPI]. As a by-product of the proof we obtain the strong Feller property of the 
associated resolvent. 
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1. Introduction. In this paper we estabhsh well-posedness of the martingale problem 
for certain degenerate second order elliptic operators. The class of operators we consider 
arises from models of catalytic branching networks including catalytic branching, mutually 
catalytic branching and hypercyclic catalytic branching systems (see [DF] for a survey of 
these systems). For example, the hypercyclic catalytic branching model is a diffusion on 
M^, d>2, solving the following system of stochastic differential equations: 

dX« = {di - X'i^)dt + ^27,(X,)xf+^^xi^)dS^ z = 1, . . . , cZ. (1.1) 
Here X{t) = {xj:^\ . . . , xj:'^^), addition of the superscripts is done cyclically so that 

Uniqueness results of this type are proved in [DPI] under Holder continuity hy- 
potheses on the coefficients. Our main result here is to show the uniqueness continues 
to hold if this is weakened to continuity. One motivation for this problem is that for 
cZ = 2, (1.1) arises in [DGHSS] as the mean field limit of the block averages of a system 
of SDE's on a hierarchical group. The system of SDEs models two types of individuals 
interacting through migration between sites and at each site through interactive branching, 
depending on the masses of the types at that particular site. The branching coefficients 7^ 
of the resulting equation for the block averages involves averaging the original branching 
coefficients at a large time (reflecting the slower time scale of the block averages) and so 
is given in terms of the equilibrium distribution of the original equation. The authors 
of [DGHSS] introduce a renormalization map which gives the branching coefficients 7j of 
the block averages in terms of the previous SDE. They wish to iterate this map to study 
higher order block averages. Continuity is preserved by this map on the interior of M'^, 
and is conjectured to be preserved at the boundary (see Conjecture 2.7 of [DGHSS]). It is 
not known whether Holder continuity is preserved (in the interior and on the boundary), 
which is why the results of [DPI] are not strong enough to carry out this program. The 
weakened hypotheses also leads to some new methods. 

The proofs in this paper are substantially simpler in the two-dimensional setting 
required for [DGHSS] (see Section 8 below) but as higher dimensional analogues of their 
results are among the "future challenges" stated there, we thought the higher- dimensional 
results worth pursuing. 

Further motivation for the study of such branching catalytic networks comes from 
[ES] where a corresponding system of ODEs was proposed as a macromolecular precursor 
to early forms of life. There also have been a number of mathematical works on mutually 
catalytic branching ((1.1) with d = 2 and 7^ constant) in spatial settings where a special 
duality argument ([M, [DP2]) allows a more detailed analysis, and even in spatial analogues 
of (1.1) for general cZ, but now with much more restricted results due in part to the lack of 
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any uniqueness result ([DFX], [FX]). See the introduction of [DPI] for more background 
material on the model. 

Earlier work in [ABBP] and [BP] show uniqueness in the martingale problem for 
the operator A^'''"'^ on C^{R'l) defined by 

Here 6^,7^ i — are continuous functions on M^|_, with bi{x) > if = 0, and 

satisfying some additional regularity or non-degeneracy condition. If bi{x) — xjqji for 
some dx d Q-matrix (qji), then such diffusions arise as limit points of rescaled systems of 
critical branching Markov chains in which (qji) governs the spatial motions of particles and 
^i{x) is the branching rate at site i in population x = {xi, . . . ,Xd)- The methods of these 
papers do not apply to systems such as (1.1) because now the branching rates 7j may be 
zero. Although we will still proceed using a Stroock-Varadhan perturbation approach, the 
process from which we are perturbing will be more involved than the independent squared 
Bessel process considered in the above references. 

We will formulate our results in terms of catalytic branching networks in which the 
catalytic reactions are given by a finite directed graph {V, S) with vertex set V = {1, . . . ,d} 
and edge set S = {ei,...,efc}. This will include (1.1) and all of the two-dimensional 
systems arising in [DGHSS]. As in [DPI] we assume throughout: 

Hypothesis 1.1. (z,z) ^ £ for alH e and each vertex is the second element of at most 
one edge. □ 

The restrictive second part of this hypothesis has been removed by Kliem [K] in 
the Holder contimious setting of [DPI]. It is of course no restriction if |V^| = 2 (as in 
[DGHSS]), and holds in the cyclic setting of (1.1). 

Vertices denote types and an edge (i, j) G S indicates that type i catalyzes the type 
j branching. Let C denote the set of vertices (catalysts) which appear as the first element 
of an edge and R denote the set of vertices that appear as the second element (reactants). 
Let c : -R ^ C be such that for j G R, Cj denotes the unique i E C such that G S, 

and for i E C, let Ri = {j : (i^j) G £}. 

Here are the hypotheses on the coefficients: 

Hypothesis 1.2. For i eV, 

7i :i?^ ^ (0,00), br-M-i^R, 
are continuous such that |fei(a;)| < c(l -|- \x\) on M^, and bi{x) > if a;j = 0. □ 
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The positivity condition on 6^1 3:^=0 is needed to ensure the solutions remain in the first 
orthant. 

li D C M*^, Cl{D) denotes the space of twice continuously differentiable bounded 
functions on D whose first and second order partial derivatives are also bounded. For 
/ e C^(R^), and with the above interpretations, the generators we study are 

Af{x) = A^^^^^fix) = J2 lj{^>c,Xjfjj{x) + lj{x)xjfjj{x) + hj{x)fj{x). 

j&R j^R j€V 

(Here and elsewhere we use fi and for the first and second partial derivatives of /.) 

Definition 1.2.5. Let O = C(]R_|_, ]R^|_), the continuous functions from ]R_|_ to M.'^. Let 
Xt{Lj) — u>{t) for wen, and let (J-'t) be the canonical right continuous filtration generated 
by X. If is a probability on M^, a probability P on O solves the martingale problem 
MP {A, p) if under P, the law of Xq is v and for all / G Cl{M.%), 

Mf{t)^f{Xt)-f{Xo)- [ Af{X,)ds 

Jo 

is a local martingale under P. 

A natural state space for our martingale problem is 

,S=|a;eR^: (0;^ + a;^) > o}. 

(i,j)ef 

The following result is Lemma 5 of [DPI] - the Holder continuity assumed there plays no 
role in the proof. 

Lemma 1.3. If P is a solution of MP{A, v), where v is a probability on M.j^, then Xf E S 
for allt>0 P-a.s. 

Here is our main result. 

Theorem 1.4. Assume Hypotheses 1.1 and 1.2 hold. Then for any probability v on S, 
there is exactly one solution to MP{A^ v). 

The cases required in Theorem 2.2 of [DGHSS] are the three possible directed graphs 
for V = {1,2}: 

(i) £: = 0; 

(ii) f = {(2,1)} or £: = {(!, 2)}; 

(iii) £: = {(!, 2), (2,1)}. 

The state space here is S' = -{(0, 0)}. In addition, [DGHSS] takes bi{x) = Oi-Xi 
for > 0. As discussed in Remark 1 of [DGHSS], weak uniqueness is trivial if either di 
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is 0, as that coordinate becomes absorbed at 0, so we may assume > 0. In this case 
Hypotheses 1.1 and 1.2 hold, and Theorem 2.2, stated in [DGHSS] (the present paper is 
cited for a proof), is immediate from Theorem 1.4 above. See Section 8 below for further 
discussion about our proof and how it simplifies in this two-dimensional setting. In fact, in 
Case (i) the result holds for any v on aUof (as again noted in Theorem 2.2 of [DGHSS]) 
by Theorem A of [BP]. 

Although it is not required in [DGHSS], it is of course natural to ask about unique- 
ness in cases (ii) and (iii) if z/ = 5(o,o)- We do have some partial results when the process 
starts at the corner (0,0), but the regularity hypotheses are stronger than in Hypothesis 
1.2 and the techniques are quite different than those used in this paper, so we do not 
pursue this here. 

Our proof of Theorem 1.4 actually proves a stronger result. We do not require that 
the 7i be continuous, but only that their oscillation not be too large. More precisely, we 
prove that there exists £o > such that if (1.2) below holds, then there is exactly one 
solution of MP{A^ v). The condition needed is 

For each i = 1, . . . ,d and each x e M^|_ there exists a neighborhood A^^; such that 

Oscjv^7i<£o, (1-2) 

where Osc a/ = sup^ / - inf^ /. 

Our proof of Theorem 1.4 is an perturbation argument, and some of our argument 
follows along the lines of [ABBP]. The operators from which we are perturbing are now 
different and the method of [ABBP] for obtaining estimates no longer applies. This 
leads to some new methodologies. 

The analytic tool we use is Cotlar's lemma. Lemma 2.13, which is also at the heart 
of the famous Tl theorem of harmonic analysis. For a simple application of how Cotlar's 
lemma can be used, see [Fe], pp. 103-104. 

We consider certain operators Tt (defined in (2.18) below) and show that 

\\Tth<c/t. (1.3) 

We require bounds on e~'^^Ttdt, and (1.3) is not sufficient to give these. This is 
where Cotlar's lemma comes in: we prove bounds on T^T* and T^*Tg, and these together 
with Cotlar's lemma yield the desired bounds on e~^*Tt dt. The use of Cotlar's lemma 
to operators arising from a decomposition of the time axis is perhaps noteworthy. In all 
other applications of Cotlar's lemma that we are aware of, the corresponding operators 
arise from a decomposition of the space variable. The bounds on TfT* and T^Tg are the 
hardest and lengthiest parts of the paper. At the heart of these bounds is an integration by 
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parts formula which refines a result used in [DPI] (see the proof of Proposition 17 there) 
and is discussed in the next section. 

In Section 3 we give a proof of Theorem 1.4. The proofs of all the hard steps, are, 
however, deferred to later sections. A brief outline of the rest of paper is given at the end 
of Section 2. 

Acknowledgment. We would like to thank Frank den Hollander and Rongfeng Sun for 
helpful conversations on their related work. 

2. Structure of the proof. 

We first reduce Theorem 1.4 to a local uniqueness result (Theorem 2.1 below). 
Many details are suppressed as this argument is a minor modification of the proof of 
Theorem 4 in [DPI]. By the localization argument in Section 6.6 of [SV] it suffices to fix 
& S and show that for some tq = ro{x^) > 0, there are coefficients which agree with 
7j, bi on B{x^, To), the open ball of radius tq centered at and for which the associated 
martingale problem has a unique solution for all initial distributions. Following [DPI], let 
Z = {i e y : x° = 0}, iVi = U^ezncR^, Ni = NiU{Zn C), and N2 = V - Ni. Note that 
iVi n Z = because x^ e S. Define 

{Xj'jjix) ifjeNi; 
x,^.7,(x) if jG (ZnC)U(7V2ni?); 
7,-(x) ifjeiVani?^ 

and note that 7^ = Jj{x^) > for all j because x^ G S. We may now write 

^'''7/(a;) ^ ^Y ^j{x)xifjj{x)^+^i{x)xifii{x) 

ieznc jeRi 

Let 5 = 5{x^) = miujg^ bi{x^) > 0, and define 

^^^'^)-\b,{x)yl -dj^N,, 

and let 6° = bj{x^), so that &° > for j ^ Ni. Although 5j(a;°) < is possible for 
J G A^2 n (and so bj may differ from bj here), a simple Girsanov argument will allow 
us to assume that bj{x^) > 5 for j E r\ (see the proof below) and so bj = bj near 
x^. With this reduction we see that by Hypothesis 1.2 and the choice of 5, bj{x) = bj{x) 
for X near x^. By changing b and 7 outside a small ball centered at x^ we may assume 
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7j > for all j, bj > for j ^ Ni, ^jihj are all bounded continuous and constant outside 
a compact set, and 

^o = X](||7-7°l|oo + ||6,-6?||oo) (2.1) 

is small. For these modified coefficients introduce 

M{x)^ 5^ Y^^j{x)xifjj{x)^^^i{x)xifii{x) 

ieznc jeRi 

+ J2 lA^)^3f3A^) + E (2-2) 
jeN2 jev 

and also define a constant coefficient operator 

i€ZnC j€Ri 

+ Y.^%fn{^)+^'f3{^) (2-3) 

i€ZnC j€N2 

As b'j < and bj\j:-=o < is possible for j G A''! (recall we have modified bj), the natural 
state space for the above generators is the larger 

5'° = S{x°) = {x eR"^ :xj>0 for all j ^ N^}. 

When modifying 7^ and bj it is easy to extend them to this larger space, still ensuring all of 
the above properties of bj and 7^ . If uq is a probability on S"*^, a solution to the martingale 
problem MP{A, fo) is a probability P on C(M+, 5"°) satisfying the obvious analogue of the 
definition given for MP(^, v). As we have Af{x) — Af{x) for x near x^, the localization 
in [SV] shows that Theorem 1.4 follows from: 

Theorem 2.1. Assume 7j : S{x'^) (0, 00), bj : S{x^) M arc hounded continuous and 
constant outside a compact set with bj > for j ^ A''!. For j < d, let 7° > 0, 6^ e M, 
65 > ifj ^ Ni, and 

Mo = max(7°, (7°)"', I^Jl) V max(6°)-^ (2.4) 

There is an £i(iVfo) > so that if sq < £i(Mo), then for any probabihty v on S{x^), there 
is a unique solution to MP{A, v). 

Proof of reduction of Theorem 1.4 to Theorem 2.1. This proceeds as in the proof 
of Theorem 4 in [DPI]. The only change is that in Theorem 2.1 we are now assuming 



7 



bj > and > for all j ^ Ni, not just bj > on {xj = 0} for j ^ Ni and > for 
j e Z n{RUC) with 6° > for other values of j ^ Ni. If > for all j e N2, then 

the proof of Theorem 4 in [DPI] in Case 1 applies without change, and so we need only 
modify the argument in Case 2 of the proof of Theorem 4 in [DPI] so that it applies if 
bj{x^) < for some j e N2. This means > by our (stronger) Hypothesis 1.2 and the 
Girsanov argument given there now allows us to locally modify bj so that bj{x^) > 0. The 
rest of the argument now goes through as before. □ 

Turning to the proof of Theorem 2.1, existence is proved as in Theorem 1.1 of 
[ABBP] -instead of the comparison argument given there, one can use Tanaka's formula 
and (2.4) to see that solutions must remain in S{x^). 

We focus on uniqueness from here on. 

The operator is the generator of a Feller branching diffusion with immigration. 
We denote its semigroup by Ql . It will be easy to give an explicit representation for the 
semigroup associated with A} (see (3.2) below). An elementary argument shows that 
the martingale problem associated with A^ is well-posed and the associated diffusion has 
semigroup 

= n n (2.5) 

ieznc jeN2 

and resolvent R\ = j e~^^Pt dt. Define a reference measure /j, on by 

IJ,{dx) = Y\ [H ^^j x\'^'^''~^dxiX JJ^ xY^' ^dxj^ n ^» n l^j- 
ieznc jeRi jeN2 ieznc jeN2 

The norm on = L^(S'°, //) is denoted by || • ||2. 

The key analytic bound we will need to carry out the Stroock-Varadhan perturba- 
tion analysis is the following: 

Proposition 2.2. There is a dense suhspace Vq C and a K{Mq) > such that 
Rx-.Vq^ Cl{S^) for allX>0 and 

[ E [ E \MR>^f)nh] + \\^^iRxf)^^\\2\ + [ E W^jW^jh] + [E II(^a/),||2 

ieznc jeRi jeN2 jev 

< K\\f\\2 for all f eVoandX> 1. (2.6) 

Here are the other two ingredients needed to complete the proof of Theorem 2.1. 

Proposition 2.3. Let P be a solution of MP{A., i') where dv = pdji for some p E L'^ with 
compact support and set Sxf = Ep ^ e~^*f{Xt) dt^ . If 

So < (2K(Mo))-^ A i48dM^)-\ (2.7) 
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then for all X> 1, 

\\Sx\\ := sup{|5,/| : ll/lh < 1} < ^ < oo. 

Proposition 2.4. Assume {P^ : x e S^} is a collection of probabilities on C(E+, 5°) such 
that: 

(i) For each x G 5°, is a solution of MP {A, 5^). 
(a) {F^,Xt) is a Borel strong Markov process. 

Then for any bounded measurable function f on and any A > 0, 



is a continuous function in x E . 

Remark 2.5. Our proof of Proposition 2.4 will also show the strong Feller property of 
the resolvent for solutions to the original MP(^, u) in Theorem 1.4-see Remark 6.2. 

Assuming Propositions 2.2-2.4 the proof of Theorem 2.1 is then standard and quite 
similar to the proof of Proposition 2.1 in Section 7 of [ABBP]. Unlike [ABBP] the state 
space here is not compact, so we present the proof for completeness. 

Proof of Theorem 2.1. Let Qfc,/c — 1,2, be solutions to MPiA,^) where u is as in 
Proposition 2.3 and define S^f = E^^J e~-^^ f{Xt) dt^ , where denotes expectation 
with respect to Qfc. Let / e C^{S^). The martingale problem shows there is a local 
martingale M-^ satisfying 

f{Xt) = f{Xo) + M^{t)+[ Af{Xs)ds. (2.8) 

Jo 

Note that for t > 0, 

Efc(sup|M/|)<2||/|U+ / Ek{\Af{S,)\)ds 



s<t 

< 



2||/||oo + c/ EJy2 Xi + l)ds<oo, 



where the finiteness follows by considering the associated SDE for X^ and using the bound- 
edness of bj. This shows that M^ is a martingale imder Q^. Let g G Vq. Multiply (2.8) by 
^g-At integrate over t, take expectations (just as in (7.3) of [ABBP]), and set / = Rxg G 
to derive 

8^x9 = J Rxgdu + S^x{{A-A^)Rxg). 
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Taking the difference of this equation when A; = 1, 2, we obtain 

\{Sl - Sl)g\ < \\S^ - SlW \\{A - A')Rxgh < \\Sl - Sl\\soK{Mo)\\gh, 

where we have used the definition of Eq (in (2.1)) and Proposition 2.2. Set £i(Mo) = 
(2K(Mo))-^ to conclude \\Sl - Sl\\ < l\\Sl - Sl\\. Proposition 2.3 imphes the above 
terms are finite for A > 1 and so we have 

^ Q > ^2.9) 

To prove uniqueness we first use Krylov selection (Theorem 12.2.4 of [SV]) to see 
that it suffices to consider Borel strong Markov processes {{Q%)xeS" ^ ^t), k = 1,2, where 
solves MP [A, dx)-, and to show that Qf = Qf for all x & (see the argument in the 
proof of Proposition 2.1 of [ABBP], but the situation here is a bit simpler as there is no 
killing). If S)^ are the resolvent operators associated with Qk, then (2.9) implies that 

j Slf{x)p{x)dii{x) = j Slf{x)p{x)dfi{x) 

for all / e L^, compactly supported p & , and A > 1. 

For / and A as above this implies S\f{x) = Slf{x) for Lebesgue a.e. x and so for all x 
by Proposition 2.4. From this one deduces Qf = Qf for all x (e.g., see Theorem VI. 3. 2 of 
[B97]). □ 



It remains to prove Propositions 2.2-2.4. Propositions 2.3 and 2.4 follow along the 
lines of Propositions 2.3 and 2.4, respectively, of [ABBP], and are proved in Sections 5 and 
6, respectively. There are some additional complications in the present setting. Most of 
the work, however, will go into the proof of Proposition 2.2 where a different approach than 
those in [ABBP] or [DPI] is followed. In [DPI] a canonical measure formula (Proposition 14 
of that work) is used to represent and bound derivatives of the semigroups P^f{x) in (2.5) 
(see Lemma 3.8 below). This approach will be refined (see, e.g.. Lemmas 3.11 and 7.1 
below) to give good estimates on the derivatives of the the actual transition densities 
using an integration by parts formula. The formula will convert spatial derivatives on the 
semigroup or density into differences involving Poisson random variables which can be used 
to represent the process with semigroup Pt from which we are perturbing. The construction 
is described in Lemma 3.4 below. The integration by parts formula underlies the proof 
of Lemma 7.1 and is explicitly stated in the simpler setting of first order derivatives in 
Proposition 8.1. 

In [ABBP] we differentiate an explicit eigenfunction expansion for the resolvent 
of a killed squared Bessel process to get an asymptotically orthogonal expansion. We 
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have less explicit information about the semigroup Pt of A and so instead use Cotlar's 
Lemma (Lemma 2.13 below), to get a different asymptotically orthogonal expansion for 
the derivatives of the resolvent Rxsee the proof of Proposition 2.2 later in this section. 

Notation 2.6. Set d= \Z (1 C\ + \N2\ = \Nf\ < d. Here | ■ | denotes cardinality. 

Convention 2.7. All constants appearing in statements of results concerning the semi- 
group Pt and its associated process may depend on d and the constants {&j,7° : j < d}, 
but, if Mq is as in (2.4), these constants will be uniformly bounded for Mq < M for any 
M > 0. 

We state an easy result on transition densities which will be proved in Section 3. 

Proposition 2.8. The semigroup {Pt,t > 0), has a jointly continuous transition density 
Pt : X ^ [0, oo), t > 0. This density, pt{x, y) is on in each variable (x or y) 

separately, and satisGes the following: 

(a) ptiy, x) = pt{x, y), where pt is the transition density associated with with parame- 
ters = 7° and 

^ \ 6° otherwise. 

In particular 

J Pt{x,y)fi{dy) = J pt{x,y)n{dx) = 1. (2.10) 

(b) If is any nth order partial differential operator in x and < n < 3, then 

sup \Dlpt{x, y)\ < C2.8r""^'--^-^^^-i n [1 + iVi/t)"^] for all y E 5°, (2.11) 
and 

sup \D^Pt{x,y)\ < C2.8r""^''~-^~^^^-i + {xj/ty/''] for all x e 5°. (2.12) 

(c) IfO<n< 3, 

sup / \D'^pt{x,y)\dii{y) < cs.gt"". (2.13) 

X J 

(d) For all bounded Borel f :S° ^R, Ptf G C^iS^), and for n < 2 and as in (b), 

D^Ptfix) = J D^pt{x,y)f{y)dM (2.14) 

and 

\\D^ Ptf Woe <C2.8t-^f Woo. (2.15) 
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Notation 2.9. Throughout will denote one of the following first or second order 
differential operators: 



D^., j < d, XiDl.^., ieZnCJeRi, or XjDl.^.J ^ iVi. 

A deeper result is the following bound which sharpens Proposition 2.8. 
Proposition 2.10. For as above and all t > 0, 

sup / \DxPt{x,y)\fi{dy) <C2.iot~^- 



sup / \Da;Pt{x,y)\ii{dx) < C2.iot ^ 
y J 



(2.16) 



(2.17) 



This is proved in Section 4 below. The case = xjD'^.^. for j e Z fl C will be the most 
delicate. 

For D as in Notation 2.9 and t > 0, define an integral operator Tt = Ti{D) by 

Ttf{x) = j Da;Pt{x,y)f{y)i^{dy), for / : ^ IR for which the integral exists. (2.18) 

By (d) above Tf is a bounded operator on L°°, but we will study these operators on 
L?{S^,li). We will use the following well known elementary lemma; see [Ba], Theorem 
IV.5.1, for example, for a proof. 

Lemma 2.11. Assume K : /S^ x 5° — > R is a measurable kernel on such that 



I 



K{;yMdy) 



< c\ and 



\K{x, ■)\lJL{dx) 



<C2. 



Then Kf(x) = J K{x,y)f{y)fj,{dy) is a bounded operator on with norm \\K\\ < ^c\Ci. 

Corollary 2.12. {a) For any f e L^fi) and t, A > 0, \\Ptfh < Wfh ^nd \\Rxfh < 
(b) Ifge Cl{S^) n L'^{n) and A^g e L'^{^i), then t Ptg is continuous in L^i/J,). 

Proof, (a) This is immediate from Lemma 2.11 and (2.10). 

(b) By {MP{A^, v)), if < s < t, then 



\\Ptg-Ps9h^ I A^Prgdr 

J s 

< I \\PrA^g\\2dr 

J 8 



<{t-s)\\A''gh. 
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We have used (a) in the last hne. 



□ 



Proposition 2.10 allows us to apply Lemma 2.11 to and conclude 

Ti is a bounded operator on LF' with norm ||Tt|| < C2.iot~^ (2.19) 

Unfortunately this is not integrable near t = and so we can not integrate this bound to 
prove Proposition 2.2. We must take advantage of some cancellation in the integral over t 
and this is where we use Cotlar's Lemma: 

Lemma 2.13 (Cotlar's Lemma). Assume {Uj : j e Z+} are bounded operators on 
L?{lJi) and {a{j) : j e Z} are non-negative real numbers such that 

V \\U*Uk\\ < a{j-kf allj,k. (2.20) 

Then 

N CO 

II J2 II - ^ •= ^• 

j=0 i=-oo 

Proof. See, e.g., Lemma XL4.1 in [T]. □ 
The subspace T>o in Proposition 2.2 will be 

-Do = {P2-jg ■.jeN,ge Cl^S"") n L\^i), X'g e L\^i)}. (2.21) 

As we can take g E C'^ with compact support, denseness of Vq in follows from Corol- 
lary 2.12(b). To see that Vq is a subspace, let P2-Hgi £ ^^o for z = 1,2 with j2 ^ Ji- If 
gi = P2-Ji-2-J2fi'i, then gi is in by Corollary 2.12 (a) and also in C^{S^) by Proposition 
2.8(d). In addition, 

M%l|2 = ||F2-.i_2-.2^Vl|2 < M"5l||2 < OO, 

where we have used Corollary 2.12(a) again. Hence P2-jigi = -P2-j2^i where gi satisfies 
the same conditions as gi. Therefore 

P2-oigi + P2-02g2 = P2-J2 {gi + ^2) e ^o- 

We show below how Cotlar's Lemma easily reduces Proposition 2.2 to the following result. 
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Proposition 2.14. There is an 77 > and C2.14 so that if is any of the operators in 
Notation 2.9, then 

\\T:TJh < C2.i4^-'-^/Vi+^/2||/||2 and 
\\TsT:fh<C2.i4S-'-^/H-'+^/^fh 

for any <t < s <2, and any bounded Borel f e (2.22) 

Assuming this result, we can now give the 

Proof of Proposition 2.2. Fix a choice of (recall Notation 2.9), let A > 1, and for 

k e define 



-fc+i 



e-^'T, ds. 



By (2.19), Uk is bounded operator on L^. Moreover if k > j then 

ds 



U. 



/.2-^ + i /.2-'=+i 

'Ukfh= / / T:Ttfdt 



<C2.142-(^/2)('=-^)||/||2. 



dsWfh 



If /c = J a similar calculation where the contributions to the integral from {s > t} and 
{t > s} are evaluated separately shows 



|[/;[/,/||2<C2.14||/||2. 



Cotlar's Lemma therefore shows that 



TV 



Uj II < v/^2(l - 2-^/^)-^ := C(r]) for all N. 



(2.23) 



Now let / = P2-^fi' £ ^0 where g is as in the definition of T>q, and for M e N set 

h = Hm = -P2-JV(i-2-*f)fi'- Then 



D^Rxf = I e-^*Pt+2-M-Nhdt 
Jo 

= exp(A2-^-^) [l)^ e-^^'Puh du 



-N-M 

M+N 00 



k+2 



e-^'^Puhdu 



= exp(A2-^-^) [ ^ Ujh + J2^~^'' / e-^('^-^)C/,_fe(Pfe/i)dw 



fc=i 
14 



In the last line the bound (2.15) allows us to differentiate through the t integral and (2.14) 
allows us to differentiate through the iJ,{dy) integral and conclude D^Puh = T^h. A change 
of variables in the above now gives 



D^Rxf = exp(A2-^-^) [ Ujh + e'^'^UoiPkh)] . 

j=0 k=l 

So (2.23) shows that 

WD.Rxfh < exp(A2-^-^)C(?7) [||/im||2 + Yl e'^'^WPkhMh 



k=l 

< exp(A2-^-^)C(r/)(l - e-'^r^hMh- (2.24) 

Corollary 2.12(b) shows that ||/tM||2 = ||-P2--^-2--~^-^5'||2 II/II2 as M — > 00. Now let 
M ^ 00 in (2.24) to conclude 

\\D.Rxf\\2<C{v){l-e-Y'\\fh, 
and the result follows. □ 

For Proposition 2.14, an easy calculation shows that for < s < 
T;T,/(x) = J K^^}{x,y)f{y)d^{y) and T,T;/(x) = j Kf}{x,y)f{y)d^{y), (2.25) 
where 

K^^Jiy.^) = / D.Ps{x,y)b^pt{x,z)dii{x), (2.26) 

Kf}{x,y) = j D,ps{x,z)DyPt{y,z)d^i{z). (2.27) 

A simple refinement of Lemma 2.11 (Lemma 3.16 proved at the end of Section 3) will show 
that (2.22) follows from 

sup 1 1 |K«(a;,y')ll<](^,2/)|rf/^(y')rf/^(^)<C2.i4«-'-^t-'+^ 

for all < t < s < 2 and i = 1, 2. (2.28) 

This calculation will reduce fairly easily to the case N2 empty and Z (iC a singleton (see 
the proof of Proposition 2.14 at the end of Section 4 below). Here there are essentially 4 
distinct choices of , making our task one of bounding 8 different 4- fold integrals involving 
first and second derivatives of the transition density pt{x, y). Fairly explicit formulae (see 
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and 



(4.7)-(4.9)) are available for all the derivatives except those involving the unique index j 
in ZnC, and as a result Proposition 2.14 is easy to prove for all derivatives but those with 
respect to j (Proposition 4.3). Even here the first order derivatives are easily handled, 
leaving = xjDxjxj- This is the reason for most of the rather long calculations in 
Section 7. In the special case (i = 2, of paramount importance to [DGHSS], one can avoid 
this case using the identity A'^Rxf = XR\f — /, as is discussed in Section 8. 

We give a brief outline of the rest of the paper. Section 3 studies the transition 
density associated with the resolvent in Proposition 2.2 for the key special case when 
Z n C is a singleton and N2 = 0. This includes the canonical measure formulae for 
these densities (Lemmas 3.4 and 3.11) and the proof of Proposition 2.8. In addition some 
important formulae for Feller branching processes with immigration, conditional on their 
value at time t, are proved (see Lemmas 3.2, 3.14 and Corollary 3.15). The section ends 
with an elementary result (Lemma 3.16) on integral operators on L^. In Section 4, the 
proofs of Propositions 2.14 and 2.10 are reduced to a series of technical bounds on the 
derivatives of the transition densities (Lemmas 4.5, 4.6 and 4.7). Most of the work here 
is in the setting of the key special case considered in Section 3, and then at the end of 
Section 4 we show how the general case of Proposition 2.14 follows fairly easily thanks to 
the product structure in (2.5). Propositions 2.3 and 2.4 are proved in Sections 5 and 6, 
respectively. Lemmas 4.5-4.7 are finally proved in Section 7, thus completing the proof of 
Theorem 2.1 (and 1.4). The key inequality in Section 7 is Lemma 7.1 which comes from 
the integration by parts identity for the dominant term (see Proposition 8.1 for a simple 
special case). In Section 8 we describe how all of this becomes considerably simpler in the 
2-dimensional setting required in [DGHSS]. 

3. The basic semigroups. Unless otherwise indicated, in this section we work with the 
generator in (2.3) where Z nC = {d} and N2 = 0. Taking d = m + 1 to help distinguish 
this special setting, this means we work with the generator 

with semigroup Pt on the state space Sm = x ]R-|. (m e N). Here we write b = 6^+1, 
7 = 7^+1 and assume 

7° > 0, 6° e M for j <m, and 7 > 0, 6 > 0. 

Our Convention 2.7 on constants therefore means: 

Convention 3.1. Constants appearing in statements of results may depend on m and 



{6°, 7? : j<m + l}. If 



Mo = Mo(7°, b') := ma^ (7^' V (7°)-^ V |6«|) V (C+i)"^ < 00, 

i<m+l 
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then these constants will be uniformly bounded for Mq < M for any fixed M > 0. 
Note that Mq > 1. 

It is easy to see that the martingale problem MP{A^, i') is well-posed for any initial 
law v on Sm- In fact, we now give an exphcit formula for Pt. Let Xt = (xf \ . . . 
be a solution to this martingale problem. By considering the associated SDE, we see that 
is a Feller branching diflFusion (with immigration) with generator 

and is independent of the driving Brownian motions of the first m coordinates. Let P^^+i 
be the law of X^'^'^^^ starting at Xm+i on C(R+,R+). By conditioning on X'^'^^ we see 
that the first m coordinates are then a time-inhomogeneous Brownian motion. Therefore 
if It = /o ^i""^'^ ds and pt{z) = {27rt)-^/^e-'' then (see (20) in [DPI]) 

/m 
/(yi, • • • , ym, X^'^^^^) Y[ p2^oj^ {yj - Xj - bp) dyj . (3.2) 

If X = {Xi, . . .,Xm+l) = {x^"^\xm+l) G Sm, let 

— — 1 

li{dx) = x^.^dx = dx^'^'>iim+i{dxm+i)- 



Recall (see, e.g., (2.2) of [BP]) that has a symmetric density qt = ql''^{x,y) {x,y > 

0) with respect to Hm+iidy), given by 



*'^(-) = (^')-''^-K^}[E;;5r(^^Tw(^^ (3.3) 

and associated semigroup Qt = Qt'^- Let 

rt{xm+i, Vm+i, dw) = P,^^, ds e dw\x'f^+^^ = ym+i) , 

or more precisely a version of this collection of probability laws which is symmetric in 
) and such that {xm+i,ym+i) n{xm+i,ym+i,dw) is a jointly continuous 
map with respect to the weak topology on the space of probability measures. The existence 
of such a version follows from Section IX. 3 of [RY]. Indeed, Corollary 4.3 of the above states 
that if 7 = 2, then 

L{X,x,y) := J expJ^^-w^fi{x,y,dw) 

^exp{(i^)(l-AcothA)}/.(^)//.(V^) ifxy>0; 
(^)'^'exp{(i^)(l-AcothA)} ifa;y = 0. 
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Here v = ^ — 1 and I,y{z) = J2m=o 7^ r(m+i/+i) ( f ) modified Bessel function 

of the first kind of index i/ > —1. The continuity and symmetry of L in {x,y) gives the 
required continuous and symmetric version of fi{xm+i,ym+i)- A scahng argument (see 
the proof of Lemma 3.2 below) gives the required version of ft for general 7 > 0. 

Now define rt{xm+i,ym+i,dw) = qt''^{xm+i,ym+i)rtixm+i,ym+i,dw), so that 

(xm+i, Um+i) i^ti^m+ij Vm+iidw) is symmetric and weakly continuous (3-5) 

and 

j j ij{ym+l,w)rt{Xm+l,ym+l,dw)lJ,m+l{dym+l) 

= Ea;^^, {'4;{xi'^^^\ It)) for all x^+i > and Borel V' : M+ ^ M+- (3.6) 

(Weakly continuous means continuity with respect to the weak topology on the space of 
probability measures.) Combine (3.6) and (3.2) to conclude that X has a transition density 
with respect to fJ,{dy) given by 

poo m 

Pt{x,y) = / TTp270u;(yj - - b°jt)rt{Xm+l,ym+l,dw) 
Jo ,=1 ' 

= pt(x(-) -y(-), wi,0,w) =p°(x(-) wi,ym+i). (3.7) 

Moreover if we set iP — (pi, . . . , b^) e and write pf {x, y) for pt{x, y), then (3.5) implies 

pf (x, y) = Pt^° {y, x) for all x,yeSm- (3.8) 

The next result is a refinement of Lemma 7(b) of [DPI]. 
Lemma 3.2. For any p > there is a cs.2{p) such that 

J w~^ft{x, y, dw) < cz.2{x + y + t)~^t~^ for all x,y>0 and t > 0. 
Proof. Assume first 7 = 2, t = 1 so that we may use (3.4) to conclude (recall = | — 1) 
= / i^Y'^ '^^fh^^ exp{|(l-AcothA)} ifx>0; 



L(A,a;,0) 
A bit of calculus shows 



exp{|(l - AcothA)} ifa; = 0. 



Acoth A > 1, aiX) := e [0, 1] for all A > 0, (3.9) 

smh A 
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with the first inequahty being strict if A > 0. The above series expansion shows that 
Iv{oLz) < a'^Iv{z) for all 2; > 0, o; e [0, 1], and so using (3.9) in the above ratio bound, we 
get 

L(A, X, y) < L(A, x, 0) for all X,x,y> 0. (3.10) 

We have 

J w~^fi{x,y,dw) — 2p J fi(x, y, [0, w^])«~^^~^(iw 

poo 

< 2pVe / L{u-^,x,y)u-^P-'^du 
Jo 

/•OO 

<2pVe L{u-^,x,0)u-^P-^du (by (3.10)) 
Jo 



00 / \ b/2 J. 

Vsinhtt-i/ 



sinhtt 

/ W ^—X 

Q \ sinh w 



< Cr 



)0/ Z rjQ 
exp{ (w coth w — l)}w^^~^dw 

exp{ — cxw^ }w'^^~ ^ dw 

/•OO 

+ / exp{-cxw - wb/2}2^'/'^w''/^+^P-^dw 



where in the last line c > and we have used (3.9), info<w<i 



infu,>i 



w coth w—1 
w 



w coth w—1 



ci > 0, and 



C2 > 0. For a; < 1 we may bound the above by (recall Convention 3.1) 



w 



2p- 



Uw+ (e-^2«;)^/2«;2p-i dw] < ci{p), 



and for a; > 1 we may, using {2we-'")^/^ < 1 for It; > 1, bound it by 

pl pOO 

Cp J exp{-cxw'^}w'^P~-^dw + J e~'''''^w'^P~-^dw < C2{p)x~^ . 

These bounds show that J w~^fi{x^ y, dw) < c(p)(l + x)~p and so by symmetry in x and 
y we get 

J w~'''fi{x, y, dw) < c{p){l + x + y)~^ for all x,y >0. 

For general 7 and t, Xg = ^X^g is as above with 7 = 2 and b = We have 
Jq Xgds = ^ Jq du, and so, using the above case. 



<c{p)t-P{t + x + y)-P. 
19 



□ 



We observe that there exist ci, C2 (recall Convention 3.1 is in force) such that 

cim''/^-^m! < r(m + 6/7) < C2m''/^-^m! for all m e N. (3.11) 
To see this, suppose first 6/7 = r > 1 and use Jensen's inequality to obtain 
r(m + r) _ f rm-i-x dx I ^ f f ^ „m-l-x 



ml 



J r(m) - \J r{m)J 



-1 r-1 

m = m 



Next suppose r e [1, 2] and again use Jensen's inequality to see that 

= Ix'-h 



r(m + r) _ f dx 



< { I x-x'^e-^'—^^y ' = (m+1)^-^ <2m^-^ 



m! ,/ r(m + 1) 

dx 
r(m + 1) 

These two inequalities imply (3.11) by using the identity T{m + r + 1) = (m + r)r(m + r) 
a finite number of times. 

Lemma 3.3. There is a C3.3 so tiiat for all t > 0, Xm+i, Vm+i > 0: 

(a) q\^\xm+uym+i) < cs.sfi""/^ + l(6/7<i/2)(^m+i A 1/^+1)1/2-6/7^-1/2]. 

(b) For alH > 0, {x,y) ^ pt{x, y) is continuous on S"^ and 
suPa;,yes^Pt(a;,y) < ca.st""'"''/^. 

(c) E,^^,(exp(-AX,('"+'^)) = (1 + A7t)-^/Texp(-a;^+iA/(l + A7t)) for all A > -(7*)"^ 

(d) If < p < 6/7 then 

E..+.((X^+'^)-^) < c,.^[j^j^^ + l)(x^+i +t)-^. 

(e) E,^^,((xi-+^))2) < C3.3(wi +^)^• 
(f) Foranyp>0, E,^^,((/o*xi"^+')ds) <C3.2(p)t-n^ + Wi)"^- 
(g) supy>o /o°°(x - y)x^/^Dlqt{x, y)dx < C3.3. 

Proof, (a) If q{x,y) = e"^-?' Em=o T^^IT^SSaT) ' ^^^^ (lt{x,y) = ht)-^/^q{x/-ft,y/-ft) 
and it suflBces to show 



9(a^,y) <C2(l + l(6/-,<i/2)(a:Ay)i/'-^/^) for all a;,y>0. 
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(3.12) 



By (3.11) and Stirling's formula we have 



— ' TniTTll t/m)' 



m=l 



(2m)! 



< c 

< c 



1 + ^1/2-6/7 (V^^g-2y5yl 

ml J 



m=l 



l + e-(v^-^)'(2^+l)i 



/2-6/7 



where in the last line we used an elementary Poisson expectation calculation (for 6/7 > 1/2 
see Lemma 3.3 of [BP] ). If 6/7 > 1/2, the above is bounded and (3.12) is immediate. 
Assume now that p= 1/2 — 6/7 > and x >y. Then the above is at most 

c(i + e-(^-^)' iV^Vvr) < c(i + Vy'iVy^ + {V^- VyF)^-^^"^^') 

<c{l + yP + yP/^)<c{l + yP). 

This proves (3.12) and hence (a). 

(b) The continuity follows easily from (3.7), the continuity of qt{-, ■), the weak continuity 
of ft{-, ■, dw), Lemma 3.2 and dominated convergence. Using Lemma 3.2 and (a) in (3.7), 
we obtain (recall Convention 3.1) 

Pt{x, y) < c{t + x^+i + y^+i)-"^/2?5''^(a;^+i, ym+i)t~'^^'^ 

< cr^/H-'Z-^t-^/' + l(,/7<l/2)(^m+l + y^+l)''/^-^/V-/2+V2-6/7 
x{Xm+lAym+l)'/'-'/''t-'/H-^/'] 

(c) This is well-known, and is easily derived from (3.3). 

(d) The expectation we need to bound equals 



/•oo 

Jo 



a;m+lVt 



= ep / v-^-P+''/'^{v + ^t)-V7e-«^-+i/(^^+7t) dv + {x^+i V t)-P 



We have used (c) in the last line. Set Cr = sup^>Q x'^e ^ and Mi = Mq . If Xm+i > t, use 



6/7 



{v + ^t)-V7e-^+i/(-+7t) < c,/^x;:i', < {cm, + l)x^+, 
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to bound the above expression for E^;^^^ by 

On the other hand if Xm+i < the above expression for Ea;^^,((X("^+i))-^')is triviaUy at 
most 

ep f'v-^-P+bh^-'^ht-'^hdv + t-P = (Z^ + At-^- 
The result follows from these two bounds. 

(e) This is standard (e.g., see Lemma 7(a) of [DPI]). 

(f) Multiply the bound in Lemma 3.2 by qt{x^y)y^/'^~^ and integrate over y. 

(g) Let h' = 7/6. Since qt{x,y) = t~^^'^ qi{x/t,y /t), a simple change of variables shows the 
quantity we need to bound is 

sup / {x- -fy)x^ dIqi {x, ^y)dx 
y>o Jo 

= sup y e"^^ / {x-y)x^'Dl(e-'=x"')dx/r(m + b') 
2/>o'fr'n Jo 



Carrying out the differentiation and resulting Gamma integrals, we see the absolute value 
of the above summation equals 



m=0 



^ e {m + l + b'){m + b') - 2m{m + b') + m(m - 1) 

b' 



y{m + b') + 2my — ym{l 



-)1 



(m>l) 



m + 6' - 1 

< ^"""^ [^'(1 + &') + y^'l - 1 + l(m>i)m(m + b'- l)-i 

e " — -7 1 1 - 6 I — 

(m + 1)! ?T 

m=l ^ ' 



m=0 



< 6'(1 + b') + 6' ^ e 



m + 6' — 1 



+ e < C3.3. 



□ 



Now let {PjJ : X > 0} denote the laws of the Feller branching process with generator 
£0/(x) = -fxf'ix). If w e C(M+, M+) let C(t^) = inf{t > : uj{t) = 0}. There is a unique 
cr-finite measure Nq on 



Cea; = {a; e C(]R+, i?+) : a;(0) = 0, C(^) > 0, uj{t) = Vt > C{i^)} 
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such that for each /i > 0, if S'* is a Poisson point process on C^x with intensity /iNq, then 

X= I has law P^; (3.13) 

see, e.g., Theorem II. 7.3 of [P] which can be projected down to the above situation by- 
considering the total mass function. Moreover for each t > we have (Theorems II.7.2(iii) 
and 11.7.3(b) of [P]) 

n^iiv : vt > 0}) = (7t)-^ (3.14) 
and so we may define a probability on Cex by 

*^ ^ No({iv : > 0)} ^ ^ 

The above references in [P] also give the well-known 

¥;{ut > x) = e-''/^\ (3.16) 

and so this together with (3.14) implies 

utdNoiiy) = 1. (3.17) 



The representation (3.13) leads to the following decomposition of X^'^'^^^ from 
Lemma 10 of [DPI]. As it is consistent with the above notation, we will use to 
denote a Feller branching diffusion (with immigration) starting at Xm+i and with generator 
given by (3.1), under the law ^Xm+i- 

Lemma 3.4. Let < p < 1. 

(a) We may assume 

= X^ + Xi, (3.18) 

where Xq is a diffusion with generator A'q starting at pXm+i, Xi is a diffusion with gen- 
erator ■yxf'^x) starting at (1 — p)xm+i > 0, and Xq,Xi are independent. In addition, we 
may assume 

Xi(t) = / utE{di^) = J2 (3-19) 

where S is independent of Xq and is a Poisson point process on Cex with intensity 
(1 — p)xm+i^o, {^jij e N} is an i.i.d. sequence with common lawF^, and Np{t) = 'E.{{i' : 
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:t5 



m+l 



vt > 0}) is a Poisson random variable (independent of the {cj}) with mean 
(h) We also have 

/ Xi{s)ds = / / Psdsl(^i,^^o)E{di^) + / / «is l(j.,=o)S(iii^) 

= E ^.W + AW, (3.20) 

3 = 1 

/ X'r+^Ms= J2rjit) + l2it), (3.21) 
where rj{t) — jQej{s)ds and l2{t) — Ii{t) + Jq XQ(s)ds. 

Remark 3.5. A double application of the decomposition in Lemma 3.4(a), first with 
general p and then p = shows we may write 

= Kit) + Y: e]{t) + Y: e]{t), (3.22) 

where Xq is as in Lemma 3.4(a) with p = 0, {ej(t), e|(t), j, k} are independent exponential 
variables with mean (l/'yt), N^lt), Np{t) are independent Poisson random variables with 
means pXm+i/lt and (1 - p)a:^+i/7t, respectively, and (Xo(t), {e]{t),e\{t)},Nl^{t),Np{t)) 
are jointly independent. The group of two sums of exponentials in (3.21) may correspond to 
Xi{t) in (3.18) and (3.19), and so we may use this as the decomposition in Lemma 3.4(a) 
with p = 0. Therefore we may take A^o to be Nq + Ap, and hence may couple these 
decompositions so that 

Ap < Ao. (3.23) 

The decomposition in Lemma 3.4 also gives a finer interpretation of the series expan- 
sion (3.2) for qt'^{x, y), as we now show. Note that the decomposition (from (3. 18), (3. 19)), 

X(-+i)(t)=X^(t)+ Y erit), 

where Aq, Np{t) and {ej{t)} satisfy the distributional assumptions in (a), uniquely deter- 
mmes the joint law of {xi'^^^\ Np{t)). This can be seen by conditioning on Np{t) = n. 
Both this result and method are used in the following. 

Lemma 3.6. Assume (j) : R+ x Z+ — > R+. If Np{t) is as in Lemma 3.4, then 
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(a) 



E.(0(4"'+^),iVo(t))) = / J]0(2/,n)(n!r(n + 6/7))-i 



n=0 



E,(</>(xi"^+^),iVi/2(t))) = / EE U 2-Xy,/c)(n!r(n + 6/7))-^ 
Proof, (a) Set a; = in (3.3) to see that XQ{t) has Lebesgue density 

that is, has a gamma distribution with parameters (&/7, 7t). It foUows from Lemma 3.4(a), 
(3.16) and the joint independence of ({ej(t)}, X'Q{t), No{t)) that, conditional on A^o(^) = 
-^(m+i) ^ gamma distribution with parameter (n + b/^,^t). This gives (a), 
(b) Apply (3.22) with p = 1/2 to see that the decomposition in (3.18) for p = 1/2 is given 
by (3.22) with X^{t) = X'^{t) + J2j=i^ ^jit) and ej{t) = e]{t). As in (a), conditional 
on (A^Q(t), A''i/2(^)) = (jj k), xj-^'^^'^ has a gamma distribution with parameters {j + k + 
6/7, 7^). A short calculation (with n = j + k) now gives (b). □ 

Notation 3.7 Let denote any nth order partial differential operator on Sm and let 

D^. denote the nth partial derivative with respect to Xj. 

If X e C(M+, E+), e Cex, G : M^. ^ K and t > 0, let 

AtG{X,u^) = AlGiX,u^) = g(^J^ Xs + ulds,Xt + i^l)-G(^J^ X^ds^Xt^ 

AjG{X,iy\iy^)=G(^J^ X^ + + vl ds.Xt + + vf) - g(^ X^ + ds,Xt + v}) 
-G{J^ Xs + vlds,Xt + vl^+G(^j^ Xsds.X^ 



--AtG{X + v\v'')-AtG{X,v''), 
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i=l 

3 



■^^ Jo Jo 



Lemma 3.8. If f : Sm — >^ M is a bounded Borel function and t > 0, then Ptf G C^{Sm) 
and for n < 3 

P"Pt/||oo<C3.8||/||oot-^ (3.24) 

Moreover if f e Ci,{Sm), then for n < 3, 

/n 

Ar(G,,,w/)(^, v\....vnY[ dnoiiv^))] , (3.25) 

where 

Gt,x(r^)f{I,X) ^ f{zi,...,Zm,X)Y[p2^oi{zj - Xj -b'^jt)dzj. (3.26) 

Proof. Proposition 14 and Remark 15 of [DPI] show that Ptf E Cl{Sm) and give (3.24) 
and (3.25) for n < 2. The proof there shows how to derive the n = 2 case from the n = 1 
case and similar reasoning, albeit with more terms to check, allows one to derive the n = 3 
case from the n — 2 case. □ 

Recall that Qt is the semigroup of X^'^'^^\ the squared Bessel diffusion with tran- 
sition density given by (3.3). 

Corollary 3.9. If : R+ ^ R is a bounded Borel function and t > 0, then Qtg G C^(R+) 
and for n <3, 

|P"Qt/||oo < ca.slblloot-". (3.27) 
Proof. Apply Lemma 3.8 to /(y) = 5f(yrn+i)- □ 
Notation 3.10 For t,S>0, a;("^) e W^, y e Sm, 1 < j <m, I,X >0, define 

m 

^t,x(m)^y{I,X) = Wp5+2^p{xi -yi + b^t)qg''^{ym+i,X). 

Lemma 3.11. (a) For each t > and y e Sm, the functions x — > pt{x, y) and x — > pt{y, x) 
are in C^{Sm), and if denotes any nth order partial differential operator in the x 
variable, then 

\D'^Ptix,y)\ + \D^ptix,y)\ < cg.nt-"-"^-''/^ for aii x,y e Sm and < n < 3. (3.28) 
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(b) sup^ / \D^pt{x,y)\ii{dy) < Cs.gt-'' for all t > and < n < 3. 

(c) For n < 3 and t > 0, y — > D'^pt{x, y) and x — > DyPt{x, y) are in Cb{Sm)- 

(d) Forn= 1,2,3, 

P n 

Proof, (a) By the Chapman-Kolmogorov equations, pt{x,y) = K x{pf/2{Xf/2jy))- The 
result for x pt{x,y) now foUows from Lemma 3.3(b) and Lemma 3.8 with f{x) = 
Pt/2{x, y). By (3.8) it follows for y pt{x, y). 

(b) For n = 1, 2, 3, TV e N and a; e Sm, let f{y) = sgn {D^ptix, y))l{\y\<N)- Then 

J \Dyt{x,y)\l(\y\<N)df^{y) = J Dyt{x,y)f{y)dfx{y) 

= \D^Ptfix)\, 

where the last line follows by dominated convergence, the uniform bound in (3.28) and the 
fact that / has compact support. An application of (3.24) implies 



/ 



Dyt{x,y)\i{\y\<N)dfi{y) <c3.8t 



and the result follows upon letting — > oo. 

(c) 



DxPt{x,y) = J pt/2{x,z)pt/2{z,y)dii{z) 

D2pt/2{x,z)pt/2{z,y)dii{z) (3.29) 



/ 



by dominated convergence and the uniform bounds in (a). The integrability of 
D^Pt/2{x, z) with respect to iJ,{dz) (from (b)) and the fact that Pt/2{zi •) £ Ci,{Sm) allow us 
to deduce the continuity of (3.29) in y from dominated convergence. Now use symmetry, 
i.e., (3.8), to complete the proof, 
(d) If y, 2; e Sm, 5 > 0, let 

m 

P'\z) = Y\_P5{Zi - yi)q5{ym+l, Zm+l)- 
1 

/ is bounded and continuous in z by Lemma 3.3(a) (with a bound depending on y, d). Let 
n < 3. The uniform bounds in (a) and integrability of f^'^ allow us to apply dominated 
convergence to differentiate through the integral and conclude 

D^^^,PtP'\x) ^ J j D-^^^p,{x,z)P'\z)d,,{z) 

^L»:^^3(a;,y)asnO. (3.30) 
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In the last line we have used (c). Now note that ii I,X > 0, then 

P m 



3-- 

m 



=Qdiym+u X)'[[ps^2^oj(yj - Xj - b^jt) = Gl^^^-,^y{I,X). (3.31) 
Use this in (3.25) to conclude that 

^L+.^t/^''(^) = Ew.(/ / KGl^,^,jX,u\...,u-)f[dno{i^')). (3.32) 

Combine (3.32) and (3.30) to derive (d). □ 

Lemma 3.12. (a) For each t > and y e M_|_, the functions x — > qt{x, y) and x — > gt(|/, x) 
are in (R+), and 



\Dyqt{x,y)\ < cz.i2t~''~^^^[l + ^/yfA for all x,y > and < n < 3. (3.33) 

(b) sup^>Q / \Dl^qt{x, y)\^rn+i{dy) < 03.12^"'^ for alH > and < n < 3. 

Proof. This is a minor modification of the proofs of Lemma 3.11 (a),(b). Use the bound 
(from Lemma 3.3(a)) 

qtM;y) < cs.sit-'^/^ + l(V7<i/2)2/'/'-''/^*-'/') 
<C3.3t-'/^(2+ V^A) 



in place of Lemma 3.3(b), and (3.27) in place of (3.24), in the above argument. Of course 
this is much easier and can also be derived by direct calculation from our series expansion 
for qt. □ 

Expectation under is denoted by and we let e(s), s > 0, denote the canonical 
excursion process under this probability. 

Lemma 3.13. If < s < t, then 

E*(e(s)|e(t)) = {s/t){{s/t)e{t) + 2-f{t - s)) < e{t) + 27s, P* - a.s. 



i 

' dx 

f e Cb(M_|_) has compact support, then Proposition 9 of [DPI] and (3.14) show that 



Proof. Let P^ denote the law of the diff'usion with generator 'yx-^ starting at h. If 



El{e{s)f{em = \imh-'El{XJ{Xt)^x,>o)ht 

= - [ [ ^"^''T'^^ Vit-siy, ^)/(^) dydx. (3.34) 
hio s J J h 
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Here we extend the notation in (3.3) by letting q^''^ denote the absolutely continuous 
part of the transition kernel for E° (it also has an atom at 0). We have also extended 
the convergence in [DPI] slightly as the functional e{s)f{e{t)) is not bounded but this 
extension is justified by a uniform integrability argument-the approximating functionals 
are bounded. By (2.4) of [BP] 

h I 7S J ■'^^ (m + 1)! V7S/ m! V7S/ 

exp(-|//7s)(7s)~-^ as I 0, 

and also ^qg''^{h,y) < (75)"^. Dominated convergence allows us to take the limit in 
(3.34) through the integral and deduce that 

E*Ms)fiem = l J J eM-yMM-'yq^,Ll{y,x)dyfix)dx. 
By (3.16) we conclude that 

El{e{sMt) = x) = {t/sfe^/^'J e-y/'^^yqtl{y,x) dy. 

By inserting the above series expansion for q^'Jgijj, x) and calculating the resulting gamma 
integrals, the result follows. 

Lemma 3.14. Let v = 6/7 — 1, and Kj^{z) = -f-{z)z + 1 for 2; > 0, where Ki/(0) = 
lim^io 'vi/(-s) = v-\-\. Then 

< C3,i4[t^ + t{Zm+l +y)]- 

Proof. Write X for X^'^'^^^ and z for z^+i- The scaling argument used in the proof of 
Lemma 3.2 allows us to assume t = 1, 7 = 2 and 8 = 26/7. Dominated convergence implies 
that 



E 



Xsds\Xi = y) = lmi^~X-'-^E,(exp(^{-Xy2) X, ds) |Xi = . 



The right side can be calculated explicitly from the two formulae in (3.4) and after 
some calculus we arrive at 
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where v = b/2 — 1 = 6/7 — 1. This gives the required equahty. 

To obtain the bound (recall Convention 3.1) it suffices to show 

K'l^iz) < c{e){l + z) for all 2; > and + 1 e [e, e~'^] 



Set a = V + 1 and recall that 



00 ^ 

/„(.) = (./2rj:--— (V2)-, 



^n!r(n + a) 

n=0 ^ ' 



and 



00 ^ 

.I'M = 2(./2)- ^-^—(.nf + 

n=l ^ ' ^ ' 

Taking ratios of the above and setting w = we see that it suffices to show 



^ (n- iy.T(n + a) 

n=l ' ' 

We claim in fact that 



— ' n\T(n + a) 

n=0 ^ ' 



for all > 0, £ < a < (3.35) 



w 



1 



(n- l)!r(n + a) " 2 



<-(a-V2vl) 



n-1/2 



+ 



n+1/2 



(n — l)!r(n — 1 + a) n!r(n + ct) 

for all n >!,«;> 0, a e [£,£"^]. (3.36) 



Assuming this, (3.35) will follow with co{s) = e"-*^/^ by summing (3.36) over n > 1. The 
proof of (3.36) is an elementary application of the quadratic formula once factors of tu""^/^ 
are canceled. □ 

Corollary 3.15. If No{t) is as in Lemma 3.4 then 



xi^+'\ No{t)) < C3.i5[t'(l + No{t)) + t{xl^^^^ + z^+i)]. 



Proof. We write X for X^'^'^^^ and z for Zm+i- Recall the decomposition (3.21): 



/ Xsds= XQ{s)ds+ / / i'sdsl(^^^^Q)E{di') + V] / ei{s)ds, 
Jo Jo J Jo j^]^ Jo 

where the second integral is independent of {{ej{t),j e N},A^o(*) = Jo ^{i/t>o)^i^^)i ^0) 
by elementary properties of Poisson point processes and the independence of Xq and S. 
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Therefore 



eJ[ Xsds 



X^(t),Aro(t),{e,(t),jeN}) 

Vo(t) t 

J] E,(/ e,(s)dsX^(t),iVo(t),{e,(t),jeN}] 



+ 



In the last hnc wc have used the independence of Xq and S, of A/o(t) and {e^, j G N}, and 
the joint independence of the {ej} (see Lemma 3.4). Now use Lemmas 3.13 and 3.14 to 
bound the last and first terms, respectively, and note the second term is bounded by the 
mean of Jq Xi(s) ds, where Xi is as in (3.19). This bounds the above by 



^^ Noit) 

C3.i4t'^ + tz + tX'oit)]+E,U Xi{s)dsj+ J] (ei(t)t + 7t^) 
< c[t^ + tXt + zt] + -fNo{t)t'^. 
Condition the above on a{NQ{t), X^) to complete the proof. □ 

We now return to the general setting of Propositions 2.2 and 2.8. 
Proof of Proposition 2.8. From (2.5) we may write 

Pt{x,y)= Yl p\{x^i),yi^i)) ql{xj,yj), 

i€ZnC j€N2 

where pi are the transition densities from Lemma 3.11 and qj are the transition densities 
from Lemma 3.12. The joint continuity and smoothness in each variable is immediate from 
these properties for each factor (from Lemma 3.11 and (3.3)). (a) is also immediate from 
(3.8). The first part of (b) is also clear from the above factorization and the upper bounds 
in Lemmas 3.11(a) and 3.12(a). The second part of (b) is then immediate from (a), (c) 
also follows from Lemmas 3.11(b) and 3.12(b) and a short calculation. 

(d) is an exercise in differentiating through the integral. As we will be doing a lot 
of this in the future we risk boring the reader by outlining the proof here and refer to this 
argument for such manipulations hereafter. Let / be a bounded Borel function on S'^ and 
< n. Note first that the right-hand side of (2.14), I{x), (finite by (c)) is continuous in x. 
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To see this, choose a unit vector Cj, set x[ = Xi ii i ^ j and Xj variable, and note that for 
h>0, 

\I(x + hej)-I{x)\< [ r \D^^Dlpt{x',y)\dx'J{y)ii{dy) 



J 

n—l 1 



< ||/||ooC2.8t-"-'/i. 

We have used (c) in the above. 

Let fniy) = f{y)'^{\y\<N)- By the integrabihty in (c), the left-hand side of (2.14) 

equals 



lim / D^pt(x,y)fNiy)dii{y) = lim / pt{x,y)fN{y)dii{y), 

N ^oo J N ^oo J 

where the differentiation through the integral over a compact set is justified by the bounds 
in (b) and dominated convergence. The bound in (c) shows this convergence is uniformly 
bounded (in x). For definiteness assume n — 2 and — D"^.^. for i ^ j. By the above 
convergence and dominated convergence we get 



^0 ^ •'V^oojO ^0 



Da:Pt{x,y)f {y)dii{y)dxidxj ^ lim / / pt{x,y)fN{y)dn{y)dxidxj 



= /[pt(a:,2/)i:!ioi:;:?]/(2/)M2/) 



Now differentiate both sides with respect Xj and then and use the Fundamental Theorem 
of Calculus and the continuity of J D'^pt{x,y)f{y)diJ,{y), noted above, to obtain (2.14). 
This shows Ptf G as continuity in x was established above. Finally (2.15) is now 
immediate from (2.13) and (2.14). □ 

We end this section with a corollary to Lemma 2.11. 
Lemma 3.16. Assume K : Sq x Sq ^ M. satisGes 



\K\ 



j j \K{x,y')\ \K{x,-Mdy'Mdx)\\^ < oo. (3.37) 



Then Kf{-) = J K{-^y)f{y)ji{dy) is bounded on L?'{ii) and its norm satisGes \\K\\ < \K\. 

Proof. Let K*K{-, ■) denote the integral kernel associated with the operator K*K. The 
hypothesis implies that 



\K*K{y,-)\^i{dy) 



< \K 



2 



By Lemma 2.11 and the fact that K*K is symmetric, we have 

IlKlP = \\K*K\\ < IkP. 
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□ 



4. Proofs of Propositions 2.14 and 2.10. 

By Lemma 3.16, Proposition 2.14 will follow from (2.28). We restate this latter 
inequality explicitly. Recall that D^; is one of the first or second order partial difi'erential 
operators listed in Notation 2.9. (2.28) then becomes 



sup 

y 



and 



sup 

y 



J J DzPsiz,x)D;,pt{z,y')iJ,{dz) iJ,{dy') J D^ps{z,x)DzPt{z,y)iJ,{dz) iJ,{dx) 

< C2.iAS~^~'^t^^'' for all < t < s < 2, (4.1) 

J j DxPs{x,z)by,pt{y\z)ii{dz) ii{dy') J D^Psix, z)DyPt{y, z)ij,{dz) ii{dx) 

< C2.i4S~^~''t^+'' for all < t < s < 2, (4.2) 



We have stated these conditions for bounded times for other potential uses; in our case we 
will verify (2.28) for all <t < s. Recall also our Convention 3.1 for constants applies to 
C2.14 and 7]. 

Until otherwise indicated, we continue to work in the setting of the last section and 
use the notation introduced there. In particular. Convention 3.1 will be in force and the 
differential operators in Notation 2.9 are 



Dx = Dxi, i <m+l, OT Da; = Xm+iD^., i <m+l. 



(4.3) 



In [DPI] a number of bounds were obtained on the derivatives of the semigroup 
Ptf', (3.24) in the last section was one such bound. Propositions 16 and 17 of [DPI] state 
there exists a C4.1 such that for all as in Notation 2.9, t > and bounded Borel function 

/, 

sup \DxPtf{x)\ < C4.it-i/||oo. (4.4) 

Although these results are stated for m = 1 in [DPI], the same argument works in m + 1 
dimensions (see, for example. Proposition 20 of [DPI]). As a simple consequence of this 
result we get: 



Lemma 4.1. For all t > 0, and all Dx as in (4.3) 



sup / \DxPt{x,y)\ii{dy) < C4^.it \ 



(4.5) 
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and 



sup [xm+i\Dl^^^qt{xm+i,y)\ + \Dx^^^qt{xm+i,y)\]l^m+i{dy) <C4.1t ^ (4.6) 

a;m+i>0 J 

Proof. Apply (4.4) and (2.14) to f{y) = sgn {DxPt{x,y)) to obtain (4.5), and to f{y) = 
sgT^{Dxm+iQtixm+i,ym+i)) to obtain (4.6). □ 

One of the ingredients we will need is a bound like (4.5) but with the integral on 
the left with respect to x instead of y. For derivatives with respect to Xi,i < m, this is 
straightforward as we now show. 

By differentiating through the integral in (3.7) we find for i <m, 

f°° (yi — Xi — Wt) -Pr 

DxiPt{x,y) = / —[lP2^^w{yj - Xj - b'-t)rt{xm+i,ym+i,dw), (4.7) 

7=1 



Dl,Pt{x,y) = 



27> 



(27»-^ 



X n^'270w(yj - ^3 - bp)rt{xm+i, ^/m+i, dw), (4.8) 



i=i 



and 



Dl,Ptix,y) = 



(j/i -Xj- h^tf 
(27» 



Zivi-Xi-h'lt) (27» 



X 



WP2'y°.w{yi - - bp)rt{Xm+l, ^m+l, dw). (4.9) 

Integration through the integral is justified by the bounds in Lemma 3.2 and dominated 
convergence. 

Lemma 4.2. For all t > 0, and i <m, 

(a) j \D^.pt{z,y)\dz^'^'> < C4.2t"^/^(t + ^m+i + ym+i)"^/^?t(ym+i, ^m+l)- 

(b) For all <p < 2, / zP^+^\D,^pt{z,y)\^idz) < C4.2t-'/'(t + l/^+i)^-'/'. 

(c) For = Zm+iDj. or D^., 



and 



sup / \D^pt{z,y)\ii{dz) < C4,2t \ 
y J 

sup J \DzPt{z,y)\n{dy) <C4.2t~^- 



(4.10) 
(4.11) 



34 



Proof, (a) By (4.7) the integral in (a) is 

■iy,-z,~b'^t) 




I n P'^i^wiVj - - b^t)rtizm+i,ym+i, dw) dz'^"''^ 

/■OO 

<c w~'^/^rt{ym+i, Zm+1, dw) 
Jo 

where in the first inequahty we used the symmetry of rt (recall (3.5)) and in the last 
inequality we have used Lemma 3.2. 

(b) Integrate the inequality in (a) to bound the integral in (b) by 

Ct J -2to+i(^ + -^m+l + J/m+l) ^^'^QtiVm+l, Zm+l)lJ'm+l{dZm+l) 

<ct-'/\t + ym+ir'/'Ey^^Aixi'-+y) 
<ct-y\t + ymr-'/'. 

In the last line we used Lemma 3.3(e). 

(c) For = (4.10) follows from (b) upon taking p = 0. The other cases are similarly 
proved, now using (4.8) for the second order derivatives. ((4.11) is also immediate from 
Lemma 4.1.) □ 

Consider first (2.28) for = Dr^. or Xm+iDl. for some i < m. 

Proposition 4.3. If = -Da;, or Xm+iDxi some i < m, then (2.28) holds with 
V = 1/2. 

Proof. Consider (4.1) for Dx = Xm+iDl We may as well take i = l. Assume <t<s 
and let 



J = 



DgPs{z,x)DzPt{z,y')iJ,{dz) iJ,{dy'). 



Then 



- J j Zm+iDl,^ps{y[,Z2,...,Zm+i,x)zm+iDl^Pt{z,y')n{dz) fi{dy') 

, 2:2, ... , Zm+i,x) - D^^psiz, x)]zm+iD\pt{z, y')ii{dz) \iJ.{dy') 



+ 

=Ji + J2. 

To evaluate Ji do the dzi integral first and use (4.8) to see 



j Zm+iD\pt{z, y') dzi = 0, 
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and so Ji = 0. (Lemma 3.2 handles integr ability issues.) 

Let J2 and J2 denote the contribution to the integral defining J2 from {zi < y'l} 
and {zi > y[}, respectively. Then by (4.8), 



J2< 




X 



^{zi<z'<y[)Dl,Ps{z' , Z2,..., Zm+l,x) 



Do the zi integral first in the above and if X = —Mi:t£^hl ^ note the absolute value of this 
integral is 



f 

J — f 



- 1 



V2^\w{y'\ -zi- bit) dzi 



= / [v^ — ^]pi{v) dv 

J —00 

<c|X|e-^'/2 

< c\y[ -z' - b1t\p2^oM -z'- b^t), 



(4.12) 



where the first inequality follows by an elementary calculation-consider |X| > 1 and \X\ < 
1 separately, note that J^{v'^ — l)pi{v) dv = 0, and in the last case use 



X 



Pi{v)[v'^ — l]dv 



\x\ 



Pi{v)[l — v'^]dv. 



Take the absolute value inside the remaining integrals, then integrate over dy2...dy'^, and 
use (4.9) to express the third order derivative. This and (4.12) lead to 



J2< 




Z. 



m+1 




l(«'<2/i) 



roo | _/_!0j3 



,0_./\-2 



X P2'y0w'{^l - Z' - b'ls)YlP2'y°w'{^j " " biS)dZjrs{Zm+l, Xm+1, dw')^ 

=2 

o P2'y2w{yi - - b'it)rt{zm+i,y'm+i,dw)dy[dz'\ 



+ ^ + 4 

'0 



X fJ'm+l{dZm+l)lJ'm+l{dy'^^i). 
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Do the trivial integral over dz2 ■ ■ ■ dzm and then consider the dy'i dz' integral of the resulting 
integrand. If z" = z' — xi + b^s and y'{ = y[ — xi — b1{t — s), this integral equals 

U"|3 




^{z"<y'^+b°t) 



+ z" 



/o /\ ' \ 11"^ J yi'yVw'K'^ J o 



X P2j°wiyi - z")dy'l dz" 

< c(27V)~'/'(27i°«;)~'/'- 
Use this in the above bound on J2 and the symmetry of rs{zm+i, Xm+i, ■) in {zm+i,Xm+i) 
(recall (3.5)), and conclude 



J'2<cj J J y ^^^'W ^^'^rs{Xm+l,Zm+l,dw')rtiZm+l,y'm+l,dw) 

X fJ'm+l{dy'fn+l)fJ'm+lidZm+l) 

-1/2 



= CJ j zl,_^:^{w') X^"^+l)dr) ys{Xm+l,Zm+udw')lIm+l{dZm+l) 

I ^'^'^~^^^^^+^^~^^'^^^ Zm+l)~^^'^rs{Xm+l, Zm+1, dw')Hm+l{dZm+l) 

(by Lemma 3.3(f)) 
< ct-V^E,^,, ((Xi-+i))'^'(£x(-+i)dr)"'^') (recall (3.6)) 

Another application of Lemma 3.3(e,f) now shows 

Symmetry (switching Z\ and y'^ in the integral defining J'^ amounts to switching the sign 
of 6°) gives the same bound on J'<1 and hence for J. This implies that the left-hand side 
of (4.1) is at most 

ct-^/^s-^/^sup y" j zl^^^\Dl^Psiz,x)\\Dl^ptiz,y)\ii{dx)n{dz) 

< ct^^/^s^^/^sup / Zm+i\Dl^pt{z,y)\ii{dz) (by Lemma 4.1) 
y J 

where (4.10) is used in the last line. 

This completes the proof of (4.1) with i] = 1/2 for = Xm+iD'^-^. The proof for 
is similar and a bit easier. Finally, very similar arguments (the powers change a bit 
in the last part of the bound on J2) will verify (4.2) for these operators. □ 

Recall the notation pt from Proposition 2.8. 
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Lemma 4.4. If and Dy are nth and kth order partial differential operators in x and 
y, respectively, then for all t > 0, k,n < 2, DyD^pt{x,y) exists, is bounded, is continuous 
in each variable separately, and equals 

J D'^pt/2{y,z)D!^pt/2{x,z)fx{dz). 

For n < 1, DyD^pt{x,y) is jointly continuous. 
Proof. From (3.29) we have 

DxPt{x,y)= / D^pt/2{x,z)pt/2{y,z)n{dz). 



Apply (2.14), with f{z) = D'^pt/2{Xj z) and 'Pt/2 in place of pt, to differentiate with respect 
to y through the integral and derive the above identity. Uniform boundedness in {x,y), 
and continuity in each variable separately follows from the boundedness in Lemma 3.11(a), 
the L-^ -boundedness in Lemma 3.11(b) and dominated convergence. If n = 1 the uniform 
boundedness of the above derivative implies continuity in y uniformly in x and hence joint 
continuity. □ 

We now turn to the verification of (2.28) and Proposition 2.10 for Dx = -Da;r„+i or 
Xm+iD'^^_^i. The argument here seems to be much harder (at least in the second order 
case) and so we will reduce its proof to three technical bounds whose proofs are deferred 
to Section 7. Not surprisingly these proofs will rely on the representations in Lemmas 3.4 
and 3.11 as well as the other explicit expressions obtained in Section 3 such as Lemmas 
3.6 and 3.14. 

Lemma 4.5. There is a C4.5 such that for all t > 0: 

(a) If -(2M2)-i < p < 1/2 then for all j <m and all y e S^, 

j zP^^,\D,^^,p,{z,yMdz) < c,.,t-'/\t + ym+ir-'/^ (4.13) 
and for all z e Sm, 

J y?n+i\Dz^+iPt{z,y)Hdy) < c^.,t-'/\t + Zm+ir-'/'. (4.14) 

(b) If0<q<2 and -{2M§)-^ <P<2, then for all j <m and all y e Sm 

J \yj-Zj\'zl+,\Dl^^Mz,yMdz) < c,.,t'^/'-\t + ym+ir+'^/'-'. (4.15) 
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(c)IfO<q< 2, then for all j <m and z e Sm, 

I \yj-zj\'^Zm+i\Dl^^^p,iz,yMdy) < c^_,[ti-^ + f^/^-^^^l^]. (4.16) 

(d) 

sup / z'Jl,\Dl^^^pt{z,y)\i,{dz) < C4.5t-'/'. (4.17) 
y J 

(e) IfO<p< 1/2, then for all (y(^), Zm+i) e Sm, 

j Zrn+iy^rn+l\Dyrn+iDl^+M^^y)\d^^"'^f'rnMdyrn+l) < C^.^t^'^ (4.18) 

and for all j < m, 

j Zm+iyl,+^\D,.Dl^^^pt{z, y)\dz^'^^fim+iidym+i) < C4.5t^-'. (4.19) 

(f) IfO<p< 3/2, then for all j < m, 

sup / zl^,\D,^Dl^^^p,{z,y)\iJi{dz) < c^.^t^-\ (4.20) 
y J 

Lemma 4.6. There is a C4.6 such that for all < t < s, 



and 



j^^ [J \D,^^^ps{z,y)\dz'^"'^jdZm+i < C4.6t/S, (4.21) 



£ [ J \Dl^^^ps{z, y)|d^("^^] dZm+1 < C4.6ts-^ < (4.22) 
Lemma 4.7. There is a C4.7 such that if 1 < p < 2, then for all t > 0, w > and 

y{m) g ^m^ 

(l(ym+l<fi'<z™+l) ^{zm+i<w<ym+'^))^m+l \D^.^Pt{z,y)\ii{dz)iim+i{dym+i) 

< C4.7[l(.<7t)t^-'+''/^ + \^y,t)t-''^w^-^'^^^hy (4.23) 




Assuming these results we now verify (2.28) and Proposition 2.10 for Dx = Dxm+i 
■i^xm+i- '^^^ analogue of Proposition 2.10 is immediate. 
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Proposition 4.8. For as in (4.3), and all t > 0, 

sup / \Da;Ptix,y)\iJ,{dy) < C4_8t~^- 

X J 

sup / \DxPt{x,y)\iJL{dx) < CA.st~^ ■ 



(4.24) 



(4.25) 



Proof. This is immediate from Lemma 4.1, Lemma 4.2(c), Lemma 4.5(a) (with p = 0) 
and Lemma 4.5(b) (with = and p=l). □ 



2^-l 



Proposition 4.9. (2.28) holds for = -Dx^+i and r] = (2M^) 

Proof. Define r] as above and fix < t < s < 2. We will verify (4.1) even with the 
absolute values taken inside all the integrals. By (4.14) with p = 0, 

In the last inequality we used r] < 1/2 (recall Mq > 1). Therefore the left side of (4.1) 
(even with absolute values inside all the integrals) is at most 



sup C4.5 

y 



D,^^,Ps{z,x)\z-l,^^{dz)t-'+^ 

D,^+,Ps{z', x)\fx{dx)] \D,^^Mz', y)Hdz'), (4.26) 



X 



sup 

y 



where we have used (4.13) with p = —rj in the last line. Now apply (4.14) to the above 
integral in x and then (4.13) to the integral in z' , both with p = 0, and conclude that 
(4.26) is at most 

as required. The derivation of (4.2) (with absolute values inside in the integral) is almost 
the same. One starts with (4.13) with p = to bound the integral in y' as above, and then 
uses (4.14) with p = —r] to bound the resulting integral in z. □ 

It remains to verify (2.28) for Dx = Xm+iD'^^^^. This is the hard part of the proof 
and we will not be able to take the absolute values inside the integrals in (2.28). 
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Lemma 4.10. For j = 1,2, J \Di^^^pt{z,y)\dz^'^^ < oo for all Zm+i > and y e 5'^, 
and 

J Di^^^Mz,y)dz^^'> = Di^^^ J pt{z,y)dz^^^ = Di^^^qt{zm+i,ym+i) 

for all Zm+i > and y G S^- (4.27) 

Proof. We give a proof as this differentiation is a bit delicate, and the result is used on 
a number of occasions. Set j = 2 as j = 1 is slightly easier. Fix y e Sm and t > 0. By 
Lemma 4.5(b,d), 

/ \Dl^+M^^y)\ + kz^+^>e)\Dl^+M^^y)\^'(dz) < oo for aU e > 0. (4.28) 



We claim 



Zm+i J Dl^^_^pt{z,y)dz^'^^ = F{zm+i) is continuous on {zm+i > 0}. (4.29) 

Note that ^(2:^^1+1) < 00 for almost all Zm+i > by (4.28). The Fundamental Theorem of 
Calculus and Lemma 3.11(a) imply that if z'^_^_i > Zm+i > 0, then 

/ 1 Dt^.Pt iz^"-^ ,4+1,2/)- Dl^^^p, (^(-) ,zm+uy) \ dz^-^^ 

-11 \Dlpt{z^'^\w,y)\dwdz^'^^ ^ as ^^^^ Zm+i or Zm+i 4+i, 

J J Zm+l 

where dominated convergence and (4.28) are used to show the convergence to 0. This 
allows us to first conclude that 



/ 



Dl^+M^^ y)M2;(-) < 00 for all Zm+i > 0, (4.30) 



and in particular, F{zm-iri) is finite for all ^m+i > 0, and also that F is continuous. 

The differentiation through the integral now proceeds as in the proof of Proposi- 
tion 2.8(d) given in Section 3 (using the Fundamental Theorem of Calculus). The last 
equality follows from (3.7) and the definition of r^. □ 

Lemma 4.11. There is a C4.11 so that for = -D^^,^^ = -D«^+i or Zm+iD^^^^, all t > 
and all y' e Sm, 

J D,pt{z,y')fi{dz)\ < C4.ii{t + y'm+i)~' (4.31) 
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Proof. Use (4.27) to see that 

J DzPt{z,y')iJ.{dz) = J D^^^^qt{zm+i,y'm+i)l^m+i{dzm+i)- 
Changing variables, we must show 

J D,qt{z,y)z'''^-Hz\ < c^.io{t + y)-\ 



(4.32) 



The arguments are the same for either choice of so let us take = for which the 
algebra is slightly easier. Let w = y/^t and x = z/^t. By differentiating the power series 
(3.3) the left side of (4.32) is then 



poo °° 

/ (7t)-^e-5] 



w 



m=0 

oo 



m!r(m + 6/7) 



(7t)-^(E^" 



w 



m=l 

CO 



m!r(m + 6/7) 

m=l ' 

<ci(7t)-ME(iV(«;)-H(^(„)>i)) + e-»], 



[mr(m - 1 + 6/7) - r(m + 6/7)]) - (7t)"^e 



1 „— 



where N{w) is a Poisson random variable with mean w, c\ satisfies Convention 3.1, and we 
have used (72 + 6/7 — !)"^ < cn~^ for all n e N. An elementary calculation (e.g. Lemma 3.3 
of [BP]) bounds the above by 



(4.32) follows. 



□ 



Proposition 4.12. (2.28) holds for = Xm+iD^^^^ and 77 = 1/2. 
Proof. Consider general < 77 < 1 for now. (2.28) will follow from 



sup 

X 



and 



sup 

X 



J I J D,ps{z,x)D,pt{z,y')i^{dz)\fx{dy') < C4.iis-^-^t-^+^ for < t < s < 2, (4.33) 




Da;Ps{x, z)Dy>pt{y' , z)ii{dz) ii{dy') < C4.iis~^~'^t-^^'^ for < t < s < 2. (4.34) 
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To see (4.1), multiply both sides of (4.33) by 

\DzPsiz, x)\ \DzPtiz, y)\ii{dz)ii{dx). 




After taking a supremum over j/, the resulting left-hand side is an upper bound for the left- 
hand side of (4.1). For the resulting right-hand side, use (4.24) to first bound the integral 
in X, uniformly in z, by C4.8S~^ and (4.25) to then bound the integral in z, uniformly in y 
by CA.%t~^. This gives (4.1) and similar reasoning derives (4.2) from (4.34). 
Next use Lemma 4.11 to see that 




Dy'Psiy', x)DzPt{z, y')n{dz) n{dy') 

< J \Dy'Ps{y',x)\c4.ii{t + y'^+i)~^i^idy') 

< [ \Dl, ps{y',x)\^i{dy') < cs-^ < 03-^-^-^+^ 



(4.35) 



In the next to last inequality we have used Lemma 4.5(b) with q = p = 0. Therefore the 
triangle inequality shows that (4.33) (with perhaps a different constant) will follow from 

J I j {DzPs{z,x) - Dy>ps{y',x))DzPt{z,y')fx{dz)\fx{dy') < 04.12^-^-^*-^+^ 

for < t < s < 2. (4.36) 

The analogous reduction for (4.34) is easier as (2.14) with / = 1 implies 

DxPs{x,y')Dy'Pt{y' , z)ii{dz) = 0. 



Use this in place of (4.35) and again apply the triangle inequality to see that (4.34) will 
follow from 



J J {DxPs{x,z) - DxPs{x,y'))Dy,pt{y',z)n{dz) n{dy') <C4,i2S ^ "^t ^+'' 



for < t < s < 2. 



(4.37) 



Having reduced our problem to (4.36) and (4.37), we consider (4.36) first and take 
r] = 1/2 for the rest of the proof. The left-hand side of (4.36) is bounded by 

j J izm+i{Dl^^^ps{z, x) - Dl,^^^ps{z^'^\y'm+i, x))zm+iDl^^^pt{z, y')i^{dz) ii{dy') 
+ J\J izm+iDl,^^^p,{z^^\y'^^„ x) - y'^^,Dl,^^^ps{y' , x)) 



X Zm+iDl^^^pt{z,y')n{dz) n{dy') 



(4.38) 
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Use the Fundamental Theorem of Calculus (recall Proposition 2.8 for the required regu- 
larity) to see that 

X Zm+i\Dl^^^pt(z,y')\ dw ii{dz)ii{dy').) 

Now recall from (3.7) that pt{z,y) = p^{z^'^^ — y^'^\ Zm+i,ym+i)- First do the 
n{dy') Hm+iidzm+i) integrals and change variables to y" = z^'^^ — in this integral to 
see that 




Ta,l < / / / / i'^{Zm+i<w<y') + ^{y'<w<Zm+i)) 



X z, 



^ +l\D^z^+rPtiy'\ Z^+l^y'm+l)\dyym+l{d^m+l)l^{dym+l) 



\Dlps{z^'^\w,x)\dz^'^Uw 
< ""^-^ P j t~^^'^w^^'^\Dlps{z^'^\w,x)\dz^'^^ w^/^-^dw 

+ 04.7^"^^ t^'''\Dlps{z^'^\w,x)\dz^'^'^dw. (4.39) 

In the last line we have used Lemma 4.7 with p = 2. Now use Lemma 4.5(d) to bound the 
first term by ct~^/^s~^/^, and use (4.22) to bound the second term by cs~^ < ct~^/^s~^/^. 
We have proved 

Note that 



Ta,2 < 




+ 



{zm+i{Dp Ps{z^'^>,y'm+i,x)-Dp ps{y',x)) 
X Zm+iDl^^^pt{z,y')ii{dz) n{dy') 
J J{zm+i - y'm+i)Dl,^^^Ps{y', x)zm+iDl^^^pt{z, y')n{dz) n{dy') 

■.=Ta,3 + Ta,4. (4.41) 

By Lemma 4.10, 

and so using Lemma 3.3(g) we have 

Ta,4: = J\ I i^rn+l - y'm+i) Z^J,liDl^^^qt{Zm+l, y'm+l)dZm+l \Dl,^^^Ps{y' , x))\fi{dy') 

< C3.3 / \Dl, ps{y',x))\fx{dy') < c^s-^ < c^r^'^s-'''^, (4.42) 
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where we have used Lemma 4.5(b) in the last hne with p = q = 0. 
For Ta 3 use the Fundamental Theorem of Calculus to write 



Ta,3 = YS^J - y'3)DjDl,^y,{y'^^^ - a;(-) + r(^M - y'M), y^i' dr] 

Now take the absolute values inside the integrals and summation, do the integral in r last, 
and for each r carry out the linear change of variables for the other (2m-dimensional) 
Lebesgue integrals: {u,w) = {z^'^'> — y'('^\y'('^) — x^"^^ + r(z*^™^ — j/'*-'"^)) (noting that 
|c^2(m)^^/(m)| < 2™|d^idw| for aU < r < 1. This shows 

Ta,3<cY^ j j J J \Uj\zl,+i\Dl^^^p^{u,Zrn+l,y'rn+l)\duf^rn+l{dZm+l)\ 

X \D^.Dl, pl{w,y'^^^,x)\ dw firn+i{dy'm+i) (4.43) 



m+1 



m 

< C 

3 = 1 



For the last inequality we have used Lemma 4.5(b) with q = 1 and p — 2. Now use Lemma 
4.5(f) with p = (for the t"^/^ term) and then with p = 3/2 (for the {ym+i)^^"^ term) to 
conclude that 

Ta,3 < c[ts-' + t-'/^-'/^] < cst-'/h-'/^ (4.44) 

Combining (4.40), (4.42) and (4.44) now gives (4.36) (with = 1/2). 
The left side of (4.37) is at most 

J J Xm+i{Dl^+^Ps{x, z) - Dl^^^ps{x, z^'^^y'^+Mm+iDl'^^^Ptiy' , z)fM{dz) ii{dy') 



+ 




Xm+i{Dl^^^Ps{x,z^'^\y'^^^) - Dl^^^ps{x,y'))y'^^^Dl.^^j)t{y',z)n{dz) fi{dy') 



■.= Tb,i+n,2. (4.45) 
The Fundamental Theorem of Calculus gives (Lemma 4.4 gives the required regularity) 




Tb,l< I I I I I {lizm+iKwKy'^^^) + ^{y'^+^<w<zm+i))^m+l 

X \Dn,Dl^^^Ps{x,z^^\w)\y'^^,\Dl,^^^pt{y',z)\ 

dw dz'^'^^ dy'^"^^ IJ,ni+l{dZm+l)fJ'm+l{dy'm+i)- 
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Re-express pt in terms of p° and set y" = y'^"^) — z^'^^ to conclude 

X y'm+llDl' Pt{y'',y'rn+l,Zrn+l))\dy''lIrn+l{dy'rn+i)l^m+l{dZrn+l) 



X Xm+i\D^Dl^^^ps {x, z^^^ ,w)\dw dz^""^ 
<C4.7 I t-'/^w'/^w'/'r-'xm+i\D^Dl^^^p,ix,z^^\w)\dz^^Uw 



+ C4.7^ / ^'''^~^^n^+l\D^Dl^^^p,{x,Z^^\w)\dz^'^^dw. (4.46) 

In the last line we used Lemma 4.7 with p=l. Now use (4.18) with p = 1/2 to bound the 
first term by C4.5C4.7t~-'^/^s~^/^. By Lemma 4.4, the second term in (4.46) is at most 



ct-^ 




'-J0 



' t'/'^\D^Ps/2iz^'^\w,z')\dz^^Uw Xm+l\Dl^^^Ps/2{x,z')\l,{dz') 



where we have used (4.21) then (4.5). (We are applying (4.21) tops/2-) We have shown 

Tb,i < dr^l'^s-^l'^. (4.47) 
For T5^25 an argument similar to that leading to (4.43) bounds Ti,^^ above by 



cxiy y y y \'^3\ym+\\'D\>^^^p\{u,y'^^^,zm+\)\du (dZm+l) 

X Xm+i\Dy,. Dl^^^p^^{w, x^+i,y'^+^)\dw /i^+i {dy'^^^) 

<cj J C4.5[l +t~^^^{y'rn+ly^'^]Xm+l\D^^jDl^^y^{w,Xrn+l,y'rn+l)^^ 

In the last line we have used the identity p^{u, y'm+i^ -^m+i) = Pt(0, v'm+i^ ^'^^ Zm+i) and 
then Lemma 4.5(c) with q = 1. Finally use (4.19) with p = and p = 1/2 to bound the 
above by c(s~^ + t~^/^s~^/^) < ct~^/^s~^/^. Use this and (4.47) in (4.45) to complete the 
proof of (4.37). □ 

Having obtained (2.28) and Proposition 2.10 for the special case A^2 null and ZdC = 
{d}, we now turn to the general case. In the rest of this section we work in the general 
setting of Propositions 2.2 and 2.14. 
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Proof of Proposition 2.14. We need to establish (2.28) (thanks to Lemma 3.16), and 
first do this for the special case when our transition density is qt = q^''^ , that is Z (1 C 
empty and N2 a singleton. Let pt be the transition density considered above with m = 1. 
Recall from (3.7) that 



J pt{x,z)dzi = J pt{x,z)dxi = qt{x2,Z2). 



(4.48) 



Let Dy^ = Dy^ or y2Dy^y^. We claim that we can differentiate through the above integrals 
and so 

j Dx^pt[x, z)dz\ = Dx2'lt{x2, Z2) for almost all 2:2 > and all x, (4.49) 
j Dx2Pt{x, z)dxi = Dx2qt{x2, Z2) for all 0:2 > and z. (4.50) 



and 



Lemma 4.10 implies (4.50). The proof of (4.49) uses / \Dx2Pt{x, z)\dzi < 00 for a.a. 2:2 > 
by Lemma 3.11(b), and then proceeds using the Fundamental Theorem of Calculus as in 
the proof of Proposition 2.8(d) in Section 3. (The stronger version of (4.49) also holds but 
this result will suffice.) 

Consider first (4.2) for g^. Let < t < s < 2. By (4.49) and (4.50), we have for all 
X2,y2 > 0, 



Dx2Qs{x2,Z2) Dy^ qt{y2,Z2) At2 {dz2 ) 



J J Dx^ps{x,z)dxi J Dy^pt{y,z)dzi H2{dz2) 



< 




Dx2Ps{x, z)Dy^pt{y, z)ii{dz) 



dxi. 



(4.51) 



Similarly, for all 0:2 > 0, 



Dx2qs{.X2, Z2)Dy'^qt{y'2, Z2)lJ,2{dZ2) IJ'2{dy2) 

Dx2Ps{x,z)dzi J Dy>^pt{y' , z)dy[ lJ,2{dz2) H2{dy2) 
Dx2Ps (x, z)by>^pt (y', z)ii{dz) ii{dy') . (4.52) 



inf 



< inf 




Integrability issues are handled by Lemma 4.10 and Proposition 4.8. The infimum in the 
second line can be omitted as the expression following does not depend on xi. Multiply 
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(4.52) and (4.51) and integrate with respect to ^2{dx2) to see that for any j/2 > 0, 




J D^^qs{x2,Z2)Dy^qt{y2,Z2)ll2{dZ2) I^2{dy2) IJ'2{dX2) 

Dx2Ps{x, z)Dy'^Pt{y' , z)iJ,{dz) / D^^p^ix, z)Dy^pt{y, z)iJ,{dz) n{dy')ii{dx) 




< 

the last by Proposition 4.12. This gives (4.2) for qt. A similar argument works for (4.1). 

Next we consider (2.28) in the general case. Write x = {{x(^i^)i^znCi (^jOjeiVa) so 
that (from (2.5)) 

Pt{x,y)= n Pt{^U)^y{j)) Yl ^ii^j^Vj)- (4-53) 

j€ZnC j€N2 

For j e {ZnC)UN2, let x-- denote x but with xy) (if j e ZnC) or xj (if j e N2) omitted, 

fij = Yli^j l^ii let Pt{xj, y^) denote the above product of transition densities but with 
the jth factor (which may be a p| or a ql) omitted. Consider (4.2) and let = -Dx(j) be 
one of the differential operators in Notation 2.9 acting on the variable j' e {j} U Rj for 
some j e Z n C. (The case / — j & N2 is considered below.) In this case (4.53) shows 
that the left-hand side of (4.2) equals 



sup 

y 




X 



H {dz^3))^j (dZj ) i^j {dy[j))i^j {dy'.) 
Dx(^)P^s{x{j),z^j))Dy^.^pl{y^j),z^j))pl{x., z.)pl{y., z-) 

H {dz{j))fij {dz-- ) fij {dx^j) )fi'- {dxj ) . 



Take the absolute values inside the two fi^{dz^) integrals (giving an upper bound) and pull 
the p^ terms out of the fij{dz(^j)) integrals. Now we can integrate the p^ integrals using 
(3.8) by first integrating over y'~., then the first z~- integral, then the x-- integral, and finally 
the second z-- integral. This shows that the left-hand side of (4.2) is at most 



sup 




X 
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In the last line we used (4.2) for (i.e., Proposition 4.12). For / = j e N2 we would use 
(4.2) for Qt, which was established above. This completes the proof of (4.2) and the proof 
for (4.1) is similar. □ 



Proof of Proposition 2.10. By Proposition 4.8 the required result holds for each pi 
factor and then using (4.49) and (4.50), one easily verifies it for qt (as was done implicitly 
in the previous proof). The general case now follows easily from the product structure 
(4.53) and a short calculation which is much simpler than that given above. □ 

5. Proof of Proposition 2.3. 

Assume P is a solution of M{A, v) where du — pdfj. is as in Proposition 2.3, and 
assume (2.7) throughout this section. Without loss of generality we may work on a prob- 
ability space carrying a d-dimensional Brownian motion B and realize P as the law of a 
solution X of 

Xi=X'o+ [ ^2UXs)XldBl+ I h{X,)ds, i^Ni, 
Jo Jo 

xi = xi + ^2jj{Xs)xidBi + J* bj{x,) ds, jeRi,ieznC. 

Set [s]n = (([ns] - l)/n) V 0, define ^j{Xs) = 7J, bj{Xs) = 6° if s < 0, and consider the 
unique solution to 

Xr = ^0 + £ ^J2^,{X[,]JX7''dBl + 

Xr^' = + ^2^,{X^,^jXT'dBl + 
Note that 

for /c > 0, on [^, ^ii] and conditional on !Fk/n, X^ has generator but with 7^, }p 
replaced with the (random) = 7(X(fe-i)/n)) = K^(fe-i)/n)- (5-2) 

We see in particular that pathwise uniqueness of X" follows from the classical Yamada- 
Watanabe theorem. An easy stochastic calculus argument, using Burkholder's inequalities 
and the boundedness of 7, 6 and Xq, shows that 

E {{X'^'y) < Cp(l + tP) for aU t > and p e N. (5.3) 

Here Cp may depend on the aforementioned bounds and is independent of n (although we 
will not need the latter). 

By making only minor modifications in the proof of Lemma 5.1 in [ABBP] we have: 
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bi{X[s]Jds, i^Ni, 

t _ 

hj{x^su)ds, j e Ri,ie znc. 



(5.1) 



Lemma 5.1. For any T > 0, sup^<y \\X^ — Xt\\ ^ in probability as n — > oo. 
For keZ+, let 

IJ,k{dx) = Y\ ^ JJ^ rfccj^a;^*^'*^' ~^c/a;j X JJ^ a;^"^'*'-' ^dxj, 
ieznc jeRi j£N2 

and let Pf{x,y) denote the (random) transition density with respect to fj^k of the diflFusion 
described in (5.2) operating on the interval [k/n, {k + l)/n]. Proposition 2.8(b) with n = 
implies that 

p'tix, y) < c^.it-"'-' n (^'^' + yT) ^ c^.it'"'' n (1 + vT)^ X, 2/ e < t < 1, 

(5.4) 

where as usual C5.1 may depend on MQ,d but not on k. We are also using (2.7) here to 
bound hi- 

Let Slf = ¥.(^ /~ e-^^fiXl") dt^ and define \\S^\\ = sup{|5^/| : ||/||2 < 1}, where 
as usual the norm refers to the fixed measure /i. 

Lemma 5.2. If (2.7) holds, then \\S^\\ < 00 for allX>0,ne N. 
Proof. It suffices to consider \S^f\ for non- negative / e LP'^IJ)- Let 



6 = sup 



(5.5) 



A bit of algebra using (2.7) shows that 



eoM^ 



2 



S < JT^P^ < 2eoMo^ (5.6) 

Mq -Eq 

Let denote expectation starting at x with respect to the law of the diffusion with 
(random) transition density p'', and let i?^ and r^{x, y) denote the corresponding resolvent 
and resolvent density with respect to jik- Then (5.2) shows that 

Slf=Y.e-^^'-¥.{E\r.U e-^'f{X,)dt)) 



/oo 
Rlf{x)p{x)d^i{x) + (E (i?^/(X,"/ J|.F(fc_i)/ J) 

k=i 

00 „ 

1 

00 

<A-^||p||2||/||2 + $]e-^'=/'^E(4"), (5.7) 
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where 



4"= / fr^xMp',J^iXl^k-i)/n,^)f^k-iidx)m 

and we have used Corollary 2.12(a) to see ||f2^/||2 < A~"'^||/||2- 

As usual we suppress dependence on d and Mq in our constants (which may change 
from line to line) but will record dependence on n. Use (5.4) and (5.5) to see that 

J r^x{^,y)p1j^{X}^k-i)/n,^)f^k-i{dx) (5.8) 
<cn f rl{x, 2/) n (1 + V^) n + ^f')^k{dx) 
<cr^ f riix^y) n {xf' + x]^^^^'),,^{dx). 

Let Ej, denote expectation with respect to the diffusion with transition kernel (x, y) — 
Pi{y,x) (as in Proposition 2.8(a)) and f\{x^y) — r^{y,x) be the associated resolvent 
density with respect to fik- Then (5.8) is bounded by 



= c. Te-^^ n E;((X:)-^^ + (X])W2.)^^, (by (2.5)) 



In the next to last line we have used the conditional independence of {X^ : i ^ Ni} (recall 
(2.5)). In the last line we have used (5.6) and (2.7) to see that 1/2 + 25 < 1, and also used 
Lemma 3.3(d). We can apply this last result because for all i ^ Ni and all A; e Z_|_, 

b^h^ - 25 > - 25 (5.9) 

Mo + eo 

> (4Mo')-^ - 4£oM3 > {8M^)-\ 

thanks to (5.6) and (2.7). For i ^ Ni we have E J(X*) <yi + Mqs, and by (2.7) and (5.6) 
we have 2d6 < 1/12. Therefore if f{s) = s''^^ + {Mos)^/'^+^^ (clearly / > 1) we may now 
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bound (5.8) by 

nOO 

nOO 

<cj e-'^f{sr'^\ dsll{l + y]/'^'') < c.,, J] (1 + vl'"^"')- 

In the definition of use Holder's inequafity and then this bound on (5.8) on one of the 
resulting squared factors to see that 



Ik < C\\J 112 



< Cn,\\\J ||2 





X ll{y!' + y-''Mdy)|,k.^{dx)\ 



<Cn,A||/||2 



X 



Here 



1/2 

U {yf' + y;''W{dyW-i{dx)] 

< Cn,x\\fh[J p1j^{Xl^k-l)/n.^)JkA^)f^k-lidx)]'^\ (5.10) 

Jk,x{^) = f rl{x,y) n {y]'"^'' + yi'')l^k{dy) 

= r ( n ((X:)V2+55 + {Xl)-'')) ds 

nE^(x:+(x:)-^^) 

e-^' n + + 

In the next to last line we have again used the independence of X*, z ^ A^i under E^, and 
the bound 1/2 + 55 < 1 which follows from (5.6) and (2.7). In the last line we have again 
used Lemma 3.3(d) whose applicability can again be checked as in (5.9). An elementary 
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oo 

As 



<c I e 



calculation on the above bound, again using (5.6) and (2.7) to see that 36d < 1/8, now 
shows that 

and therefore by (5.10), 

<Cn,A||/||2n(l + ^5-l)/n+v)- 

Take expectations in the above and use some elementary inequalities to conclude that 

<Cn,x\\fhJ2^((^+^i>^l)/n)') 

<Cn,x\\fUl + ik/n)'), 
the last by (5.3). Put this bound into (5.7) to complete the proof. □ 

Remark 5.3. The above argument is considerably longer than its counterpart (Lemma 
5.3) in [ABBP]. This is in part due to the non-compact state space in the above leading 
to the unboundedness of the densities on this state space (recall the bound (5.4)). More 
significantly, it is also because the argument in [ABBP] is incomplete. In (5.4) of [ABBP] 
the norm on / actually depends on k and is not the norm on the canonical space. The 
argument above, however, will also give a correct proof of Lemma 5.3 of [ABBP]-in fact 
the compact state space there leads to considerable simplification. 

Proof of Proposition 2.3. This now proceeds by making only minor changes in the 
proof of Proposition 2.3 in Section 5 of [ABBP]. One uses the above Lemmas 5.1, 5.2 and 
Proposition 2.2. We only point out the (trivial) changes required. For / e C^{S^) one uses 
Ito's Lemma and (5.1) to obtain the semimartingale decomposition of /(X"). The local 
martingale part is a martingale as in the proof of Theorem 2.1 in Section 2 (use (5.3)). 
Corollary 2.12(a) is used, instead of the eigenfunction expansion in [ABBP], to conclude 
that the constant coefficient resolvent Rx has bound A"-*^ as an operator on L^. The rest 
of the proof proceeds as in [ABBP] where the bound Eq < {2K{Mo))~^ is used to get the 
final bound, first on 'S'^(|/|), and then on 5^(1/1) by Fatou's lemma. □ 

6. Proof of Proposition 2.4. 
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Let {¥^,Xt) {x e S^) be as in the statement of Proposition 2.4. Throughout this 
section, for any Borel set A we let Ta = Ta{X) = inf{t : Xf e A} and ta = TAiX) = 
mi{t : Xt ^ A}, he the first entrance and exit times, respectively, and let \A\ denote the 
Lebesgue measure of A. We say a function h is harmonic m D = B{x,r)r\S if h is bounded 
on D and h{Xt/\T-^) is a right continuous martingale with respect to Fx for each x. 

The key step in the proof of Proposition 2.4 is the following. 

Proposition 6.1. Let z E S. There exist positive constants r, ce.i and a, depending on 
z, such that if h is harmonic in B{z, r) fl S, then 

\h{x)-h{z)\ <C6if ^^"^^ ^"f sup \h\\, X e B{z,r/2)nS. (6.1) 

\ r / ^B{z,r)ns ' 



Proof. By relabeling the axes we may assume that = {a; e : a^j > for z > Jq}. 
If z is in the interior of 5°, the result is easy, because the generator is locally uniformly 
elliptic, and follows by the first paragraph of the proof of Theorem 6.4 of [ABBP]. So 
suppose z G dS^"^ . Then Jq < d and we may assume, again by reordering the axes, that 
there is a, K E {Jq, . . . , d—1} so that Zi = for all i > Kand Zi > if Jq < i < K . Assume 
(set min0 = 1) 

< r < min — . (6.2) 

Jo<i<K 2 

Since our result only depends on the values of of h in B{x, r) fl we may change the 
diffusion and drift coefficients of the generator of X outside B{z, r) D S as we please. By 
changing the coefficients in this way and again relabeling the axes if necessary, we may 
suppose that our generator is 

J K 

(Ji (x)Xq( 

i=l i=J+l 

d d 

+ Y ai{x)xifii{x) + ^bi{x)fi{x), (6.3) 

i=K+l i=l 

where J < K, a{i) e + . . . , for z = J+1, . . . , -R', each cTj is continuous and bounded 
above and below by positive constants, each bi is continuous and bounded, and each hi 
foTi>K is bounded below by a positive constant. We have extended our coefficients to 
the possibly larger space = {x E : Xi > for all i > K} as this is the natural state 
space for A. As B{z,r) f] = B{z,r) fl (by (6.2)) this will not affect the harmonic 
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functions we are dealing with. For < 5 < 1 let 

J K d 

QniS) = l[[zi - 2-/^ z, + 2-"/2] X n - 2""' + 2"1 X n 2'^] 

i=l i=K+l 
J K 

Rn{5) = Y[[zi - 1 • 2-/^ + 1 • 2-"/2] X Y[[zi-i 2-^ + § • 2-] 

X JJ [(52-", 2-"]. 

Take n > 1 large enough so that (5n(0) C B{z, r/2) n 5^ 

We will first show there exist C2, (5 > independent of n such that 

^x{TR^^,i5) < TQ„(0)) > C2, X e gn+l(0). (6.4) 

We may assume there exist independent one-dimensional Brownian motions Bl such that 

where dA\ = hi{Xt) dt and 

dMi = {2a,{Xt)y/^ dBl i<J, 

dMi = (2a,(Xt)X"^'V^^ t^^t, J + 1 < i < 
dMi = {2ai{Xt)Xiy/^ dBl K+l<i<d. 

Since the bi and ctj are bounded, there exists to small such that for all x e Qn+i(0) 

P,( sup |Mj|>l-2-/2)<^, z<J, 

s<to2-" 

P,( sup |M:|>i-2-'^)<^, J + l<i<d, 

s<to2-" 



and 



sup |A*|<i-2-", l<z<d. 

s<to2- 



The first and last bounds are trivial, and the second inequality is easily proved by first 
noting 

sup / X*ds)<c2-2", 

and then using the Dubins-Schwarz Theorem and Markov's inequality. Hence 

sup P,(TH„+,(o)<to2-")<i (6.5) 
a;eQ„+i(0) 

55 



By Lemma 6.2 of [ABBP] there exists 6 such that if U is uniformly distributed on 
[to/2, to], then 

sup F{Xij <S)< 1/M, K+l<i<d. 

Scahng shows that 

sup P^(X^2-" < <^2~'') < l/4d, K+l<i<d. 

Therefore by (6.5), for any x e Qn+i(0) with P^-probabihty at least 1/2, X[/2-" £ -Rn+i(<^)- 
Since i?n+i(0) C Qn(0) and C/2-" < to2~'', this and (6.5) proves (6.4). 

Take 5 even smaller if necessary so that \QniO) — QniS)\ < ^\QniO)\ and 5 < ai < 
\bi\ < 5-^ for all i. Next we show that if G C Qn(0) and \G\ > \Qn{0)\/3, then there 
is a 03(6) > 0, independent of n, so that 

Pcc(Tg<tq„(5/2))>C3, a;ei?n+i(5). (6.6) 

Let 

yi _ I 2"'-^2-"t' i > J, 

* \ 2"/2x^_„^, i < J. 

It is straightforward (cf. [ABBP], Proof of Theorem 6.4) to see that for t < tq^(^sj2)i 
Yo e Qo(5/2), Yt solves 

dY^ = d,(Yt)dBl + h(Yt)dt, 

where the are independent one-dimensional Brownian motions, the hi are bounded 
above, and the ai arc bounded above and below by positive constants depending on 5 but 
not n. (6.6) now follows by Proposition 6.1 of [ABBP] (with a minor change to account 
for the fact that \G fl Qn{5 /2)\/\Qn{5 /2)\ is greater than rather than i). 

Combining (6.4) and (6.6) and using the strong Markov property, we see that if 
C4 = C2C3, then 

Px(Tg < rQ„(o)) > C4 > 0, x^Qn+M- 

Suppose h is harmonic on Qj^^ for some uq- Our conclusion will follow by setting a 
= log(l/p)/log2 if we show there exists p < 1 such that 

OscQ^^^(0)/i <pOscQ^(0)/i, n>no, (6.7) 

where Osc Ah = sup^ h — inf a h. Take n > uq and by looking at 05/1 + cg, we may suppose 
supg |-g-) h — 1 and iufg^^o) h = 0. By looking at 1 — /i if necessary, we may suppose 
1^1 > l\Qn{0)\, where G ^ {x e Qn{0) ■ h{x) > 1/2}. By Doob's optional stopping 
theorem 

h{x) > E^[h{XTay,TG < TQ„(o)] > iPx(TG < TQ„(o)) > C4/2, X G Qn+l(0). 
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Hence Osc q„_,_i(o) h < 1 — C4/2, and (6.7) follows with p = 1 — C4/2. 



□ 



Proof of Proposition 2.4. We can now proceed as in the proof of Theorem 6.4 of 
[ABBP]. To obtain the analogue of (6.14) in [ABBP], we note from (2.2) that if a; e dS^, 
at least one coordinate can bounded below by a squared Bessel process with positive drift 
starting at zero. 

Remark 6.2. Essentially the same argument shows that if for each x E S, is solution 
of MP{A,5x) as in Theorem 1.4 (it will be Borel strong Markov by Theorem 1.4), then 
the resolvent S\ maps bounded Borel measurable functions to continuous functions. After 
localizing the problem, one is left with a generator in the same form as (6.3) and so the 
proof proceeds as above. □ 



7. Proofs of Lemmas 4.5, 4.6 and 4.7. 

We work in the setting and with the notation from Sections 3 and 4. Recall, in 
particular, the Poisson random variables Np{t) from Lemma 3.4. 

Lemma 7.1. There is a C7.1 such that for all < q < 2, 1 < j < m, y e Sm, < t, and 

z' = Zm+i/lt > 0; 

(a) If for X > and n e Z_|_, 

M^',n) = (/ + 1)^/2-1 

and 



l(n<l) + l(n=l)^' ^ + l(n>2)2;' ^ + {j^ ' ^ 



M^', ^, n) =l(n<2)(l + ^'"")(1 + Z"^' + X'^/') 

+ lin>3)z'~\\n - z'\^ + 3n\n - z'\ + n){z"''^ + n^/^ + x'^'^), 



then 



JlVj - Zj\'^\Dl^^^Pt{z,y)\d^ 



(m) 



(7.1) 



6^0 

(b) If for X, n as in (a). 



0(z', X, n) = l(„<i)(l + X'^/' + {zr^') + ^r^=l){zT\^ + (zT^' + x'^') 

+ ^n>2){zT\ri +{n- z'f){{zr/' + n«/2 + W^), 
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then 



J \yj - Zj\'^\Dl^+^Pt{z,y)\dz 



(m) 



o— >0 



(c) 



j \Dzm+iPt{z,y)\dz 

< cr.it-' liminf E,^^, ((1 + z')-' + ^^^^)) • (7-3) 

In addition for all z e Sm, 

J\D,^^Mz,yW^^ 

<cr.it-H\^ME,^^^(^qs{yrn+i,xt^'^)(^^^ (7.4) 

Proof. The proof of (a) is lengthy and the reader may want to first take a look at the 
simpler proof of (c) given in Section 8. 

(a) By Lemma 3.11(d), Fatou's lemma and symmetry we have 

< liminf I \yj - Zj\'^\e,^^, (^A^,GI,^,JX,u\u^u^)[1^^,^^ fo, .^^^, 3^ 
+ 31(i^i>o,i/f =i.3=o) + 31(^i>o,i.2>o,,.3=o) + l(i.i>o for i=i,2,3)] 

3 

xJjNo(rfi^*)) dz^"^^ 

i+l 

:= hm inf Ef + 3E^ + 3E| + eI (7.5) 

(5^0 

Consider Ef first. An explicit differentiation shows 

\Dtpu{w)\ < cp2u{w)u-'' for A; = 1,2,3, (7.6) 

which implies 

m 

\D'iGl^,^,JI,X)\<C2qs{ym+i,X)r^l[p,s+4^oj{z,-y, + b''^t), for A; =1,2, 3. (7.7) 
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Use this (with k = 3) with the Fundamental Theorem of Calculus to see that on {ul > 
for z = 1,2,3}, 

m 3 

<crt.(y^+i,X^+'))(t^+ + + fi^'sdsY') n f^^'ds. (7.8) 
(3.14) and (3.16) imply 

j lyP No{du) = {-fsf-^Tip + 1) for p>0, (7.9) 

and so by Jensen 

= r(2 + q/2)t''/^ I {-fsy/^ds < ct''+\ (7.10) 
Jo 

This bound, (7.9) with p = 1, and (7.8), together with the expression for E(, shows that 

E( < cE,^^,(/-3g,(y„+i,X^+^))(t'^ + 5^/2 + /,''/^)t3) 

<ct^-'E,^^,(g,(2/^+i,X;-+^)))(l + (5V2/^).)(l + ^')./2-3^ ^^^l) 

where Lemma 3.2 is used in the next to last inequality. 

Let us jump ahead to which will be the dominant (and most interesting) term. 
We use the decomposition and notation from Lemma 3.4 with p = 0. Let Sn = Yl7=i 
Rn = ElLi ^iit)^ Pk{z') = e-''^ = P{No{t) = k), and iV^ = N{N -1) . . .{N -k + 1). 
What follows is an integration by parts formula on function space. Recalling that No(j>'t > 
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0) = (7^) ^ from (3.14), we have from Lemma 3.4 and the exponential law of //^ under P^* 
(recall (3.16)) that 



E 



-3 



^tGl^-^jX, v\ v\ v"-) n(l(.i>o)No(ci^^))^ 

1^2^+1 (G't,2(m)^y(-RjVo+3 + -f2(*), 'S'aTo+s + ^o(*)) 

- 3G'^^^(^)^^(i?ArQ+2 + Sno+2 + X'Q{t)) 

+ 3G^^^(m) ^y{RNo+i + hit)-, Sno+1 + ^oit)) 

^ ( 

oo 

= (7i)"'| J](Pn-3(-2') - Spn-2iz') + 3pn-liz') - Pniz')) 



n=0 



X E,^^, (Gj,,(^),,(i?n + /2(t), 5n + X^(t))) 

00 

{-ft)-^\j2Pn{z'){z')-'^[n'^^^ - 3n^2)/ + 3n(/)2 - {z'f 



n=0 



^-3 



E 



2m+l 



([ivf - 3ivi^).' + 3iVo(.')^ - izr]GUr.,,y{Iu Xt^'^))\. (7.12) 



In the last line we have again used Lemma 3.4 to reconstruct {If, xj:^~^^^). 
We also have 

J\y,-z,\'^Gl,^,Jh,xt+'^)dz(-^^ 

= qs{ym+i,X^'^^^^) J \yj - 2j|%+270/,(% - Vj + h^jt)dzj 
< cqsiVm+u X^+'^)[t^ + 5'^/' + (7.13) 
Combine (7.12) and (7.13) to derive 

Ei < cz-l,E,^^, (|ivf ) - 3<)/ + 3No{zr - {z'^l 

X q5{ym+i,xl-+'^)[t'^ + 5'^/' + if']). (7.14) 

Apply Jensen's inequality (as q/2 < 1) in Corollary 3.15 to see that 

E,^^Alf"\No,xi'^+'^) < c[t^ + tmf' + t^/\X^'^+'Y' + t'^'^d]^ (7.15) 
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and so 



X [1 + (s'/^ty + ixi"'+'^ /ty/^ + (zy/^ + (7.16) 

Next consider E2. Not surprisingly the argument is a combination of the ideas used 
to bound Ef and E^. Define 

^' =^t(-),. + + ds, X) - G^,^,^^ (/ + J* ds, X) 

- GU^\y {} + r X) + G',^,,(.),,(/, X). 

«/ 

Now apply the decomposition in Lemma 3.4 with p = 1/2 so that -/Vi/2(^) is Poisson 
with mean z' j2. Arguing as in the derivation of (7.12), but now with a simpler first order 
summation by parts (which we leave for the reader), we obtain 



(/ A?G^^^^(^)^^(x,^.^^.^^.3)l(,l>o,.-.-o) n^o(^^')) 



E 




HU-^,y ( / + -Ids, + .ly, .3) 

- ^t(^),. ( r Xi-^+'^ds, Xl^^'\ P*(ci^^)lK-.3=o) n ^oidv^)) 

•'^ k=2 

-^M(-),y(-^2(^) + Rni/2 + 1i ^^{t) + S'aTi/s + I, i^^, l^^) 

3 

((iVl/,-^72)i^f^^(^,^^(/(^),x("^+^^^.^z.3)) J[No(di^^)|. (7.17) 




k=2 



Now use (7.7) (with k = 2) and argue as in (7.8) to see that 



\yj-z^ml,^,JI,X,i.^i.')\dz(-^^l[No{du'^) 




k=2 
3 ,t 



3 



< ci-\s{ym+i,x)[t'i + 5^/2 + 
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(7.18) 



where the last hne uses (7.10). 

Take the absolute values inside the inside the integral in (7.17), multiply by \yj—Zj\'^, 
integrate with respect to z^'^^ , and use the above bound to conclude that 



X 



(t^ + 5'^/2)/(t)-2+/(t)9/2- 



(7.19) 



If r > 0, the independence of Xq from (A^i/2, {^j}) and Lemma 3.2, applied to Xq, imply 
that 



IE..+.(/(t)-nX^+^\iVv2) 



<E E / Xl,{s)ds] |X^(0,iVi/2,{e,} 



X 



(m+l) 



iVi/2) 



EE 



X 



(m+l) 



1/2. 



X',{s)ds) X'.it) 

<c{t + Zra+l/2)-H-\ 

The last line is where it is convenient that p = 1/2 > 0. 

Use (7.20) in (7.19) with r = 2 and 2 - q/2. After a bit of algebra this leads to 

Ei < ct^-3(^')-'lE,^^,(|iVi/2 - z'/2\qs{ym+r,x'r^'^)) 

, 1 



(7.20) 



The argument for is similar to the above. One works with 



(7.21) 



H: 



(/, X, u') = Gl^,^,^^ (/ + f u^ds, X) - Gl^,^,JI, X). 

J 



The required third order difference of G^^(^) ^ on {y\ > 0, > 0, = 0} is now a second 



order difference of _ _ . Minor modifications of the derivation of (7.21) lead to 



-3/'^/N-2t 



< ct'^-\zY'^,^^,{\Nfl>, - Ny^z' + {z'/2f\qs{ym+i,xl^^'^)) 

X [(^A)''(l + -s')"' + (1 + (7.22) 
The above bounds in Ef i = 1,...4 may be used in (7.5) and after the terms 
involving VS/t are neglected (for g = these terms are bounded by their neighbours, and 
for > 0, if they do not approach 0, the right side below must be infinite) we find 

J\yj-z,nDl^^Mz^y)\dz^^^ 

< ct«-3liminf E,^^, (qs{ym+uXi^+'^) \{1 + ^^1^-^ (7.23) 
X ((1 + z')-'' + |iVi/2 - zl2\{\ + z')-\zY^ + liV^^J) - Ny^z' + {z' I2f\{z')-'') 

+ \Nf^ - mf^z' + m^i^z'f - {^z')%z'Y\i + (xi"^+'V*)'/' + (^')'/' + ^^0^']])- 
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The required bound follows from the above by a bit of algebra but as the reader 
may be fatigued at this point we point out the way. Trivial considerations show it suffices 
to show the following inequalities for no,ni/2 £ Z+ and z' >0: 



\n 







3n'^''z'-3no{z'y-{z'r\{z'y 



< c[l(n„<2)(l + {z')-"') + l(no>3)(-s')" (1^0 - z'\^ + 3no|no - z'\ + no), (7.24) 



and 



[(1 + z'r' + \ni/2 - z'/2\{l + z')-\z')-' + \n[% - n./^z' + {z' /2f\{z')-^ 

z' 

^ C[l(ni/2<1) + l(ni/2 = l)^'~^ + ^im/2>2)z'~^{ni/2 + (^1/2 - -)% (7.25) 

(7.24) is easy. (7.25) reduces fairly directly to showing that for ni/2 > 2, 

(1 + z'r' < c{n^l2 + (ni/2 - z' l2f){z'Y'' ■ 

If < 1 this is trivial and for z' > \ consider ni/2 < z' jA, and ni/2 > z' jA^ separately. 
This completes the proof of (a). 

(b) The proof of this second order version of (a) is very similar to, but simpler than 
that of (a). One now only has a second order difference and three E'f terms to consider. 
In fact we will not actually need g > in (a) but included it so that the reader will not 
complain about the missing details in the proof of (b) (where g > has been used in 
Proposition 4.12). We do comment on the lack of A^i/2 in this bound. 

An argument similar to that leading to (7.23) shows that 
/ \zi - yj\'^\Dl^+iPt{z,y)\dz^'^^ is bounded by 

cf^-' liminf E {qsiVm+u ^^+'^) [(1 + zT'^'^/' ((1 + zT' + |iVi/2 - z' /2\{zr') 
+ |Arf ^ - 2Noz' + {z'f\{z')-\l + (X^^+'V*)'/' + (zT^^ + K^^)]) 
= ct«-2 hminf E (qsiVm+i, X^r^'^)[Ty2 + Tq]) . 

It is easy to check that 

To<ccl>{z',xi^+'^/t,No), 
and, using N1/2 < Nq from (3.23), that 

Ti/2 < 2(1 + z')-'+'^/'(l + Noiz'r' + 1) = ^1/2. 

Hence to prove (b), it remains to verify 

T,/2<ccf>{z',xi"'+'^/t,No). 

63 



Trivial considerations reduce this to showing that (1 + z'^-'^^il'^ < c4>{z',X^'^ '/t,No) 
for Nq > 2. This is easily verified by considering A^o < -^'/^ and A^o > -^'/^ separately. 

(c) Note that (7.3) is the first order version of (a) and (b) with q = 0. The proof 
is substantially simpler, but, as it plays the pivotal role in the proof for the important 
2-dimensional case, we give the proof in Section 8. (7.4) then follows immediately since 
the spatial homogeneity in the first m variables, (3.7), implies 

Pt{z,y) ^pt{-y^'^\zm+u-z^'^\ym+i). (7.26) 

□ 



Proof of Lemma 4.5(b). Let J denote the integral to be bounded in the statement of 
(b), and Pn{w) = e~'^w'^/n\ be the Poisson probabilities. Let r„ be a Gamma random 
variable with density 

gn{x) = a;'^+^/^-^e-"r(n + b/^)-\ (7.27) 

and recall z' = Zm+i/lt. By integrating the bound from Lemma 7.1(b) in Zm+i (using 
Fatou's Lemma) we see that 

J < cr.it^-^liminf J zf,^,E (g,(y^+i, X^+^))0(^', xi"^+^V^, ^Vo)) z'^l-'dzm+i- 

Our formula for the joint distribution of {xj:^^^\ Nq) (Lemma 3.6(a)) allows us to evaluate 
the above and after changing variables and the order of integration we see that if y' — y/^t, 
then 

/oo 
{q6{ym+i,y) ^Pn{y') 
n=0 

POO 

+ l(,>2)(z')"'[(^ - ^'f + n][iyr^' + (zr^^ + n'^l\z'Ydz^y'^l^-Uy 

/oo 
Qs{ym+i,y)y^^Pn{y') 
n=0 

X 3(n<i)E ((1 + (y')^/' + n/')ri;) + l(n=i)E ((1 + (y')«/' + ^'rl'')^^^^'') 
+ l(n>2)E (((n - rj^ + n)r-2+^'((y')^/' + r^/' + n^/^))]^''/^-^^. (7.28) 
There is a constant cq (as in Convention 3.1) so that 

_3 

E (r;) < co(n V 1)'^ for aU |r| < 4 and n e Z+ satisfying r + n > — 3 . (7.29) 



4Mo^ 
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Indeed the above expectation is T{n + 6/7 + r)/T{n + 6/7), where 

-3 



r + n + 6/7 > ^ + Mq-^ = (2Mo)-^ 

The result now follows by an elementary, and easily proved, property of the Gamma 
function. 

Assume now the slightly stronger condition 

Irl < 3, n G Z_|_ and r + n > 5-. (7.30) 

Then r„ = Fq + S*^, where is a sum of n i.i.d. mean one exponential random variables. 
If s and s' are Holder dual exponents, where s is taken close enough to 1 so that the 
conditions of (7.29) remain valid with rs in place of r, then 

E ((r„ - nfVl) < E ((r„ - nf^'y/^'E (T^) V« 

<cn(nVl)^ (7.31) 

where we have used an elementary martingale estimate for \Sn — nl and (7.29). Here c 
again is as in Convention 3.1. 

We now use (7.31) and (7.29) to bound the Gamma expectations in (7.28). It is 
easy to check that our bounds on p and q imply the powers we will be bounding satisfy 
(7.30). This leads to 

/oo 
qs{ym+i,y) ^Pn{y') 
n+O 

X "l(,<i)(l + ivr/') + \n>2)n-'+niyr^^ + n^/^)y/^-' dy 

<ct^-2+P liminf / qs{ym+i,y)\e-y'{l + y'){l + {yr/') 

+ E (l(iv(,')>2) {{y'r^'Niy'r-' + Niyy/'-'+P))] y^'/^-'dy. (7.32) 

In the last line N{y') is a Poisson random variable with mean y'. Well-known properties 
of the Poisson distribution show that for a universal constant C2 

E (Ar(yO"l(Ar(y')>2)) < ^Ay') = 02(1 + y'Y for aU y' > 0, |r| < 2. (7.33) 

For negative values of r see Lemma 4.3(a) of [BP] where the constant depends on r but 
the argument there easily shows for r bounded one gets a uniform constant. If 

h{y') = e-y'{l + y'){l + (y')^/') + K/2-i+,{y') + y"^/%-i{y% 
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then clearly 

As all of the powers appearing in (7.32) satisfy the bounds in (7.33), we may use (7.33) to 
bound the left-hand side of (7.32) and arrive at 

J < ct^-2+^'liminfE,^^,(/.(x(-+^V7t)) 

o — >U 

< ct^-2+niminf E,^^,((l + {Xf-^'^ hty'^-'^n 

= ct'^-2+P(l + ^')9/2-i+P (Dominated Convergence) 

= ct'^'^-\t + y^^,Y/''-'+^. 

□ 

Proof of Lemma 4.5(c). The spatial homogeneity (7.26) shows the integral being 
bounded equals 

j IVj - (-2;j)|^2;m+l|^L+i^'t(2/^'"^^rn+l,-^^'"\ym+l)|c^y^'"^/im+l(c^ym+l)^ 

This shows we can again use the upper bound in Lemma 7.1(b) to bound the integral over 

y{m) in 

the above. One then must integrate the resulting bound in ym+i instead of Zm+i- 
This actually greatly simplifies the calculation just given as one can integrate ym+i at the 
beginning and hence the qs term conveniently disappears (see the proof of (4.14) below). 
For example, if wc neglect the insignificant n < 1 contribution to (J) in Lemma 7.1, the 
resulting integral is bounded by 

cf^-'iz')-'^ (i(Aro>2)(iVo + (ivo - z'rwr^' + K'^ + (X(-+^V^)'/'))- 

This can be bounded using elementary estimates of the Poisson and Holder's inequality, 
the latter being much simpler than invoking Lemma 3.6. We omit the details. □ 

Proof of Lemma 4.5(a). (4.13) is the first order version of Lemma 4.5 (b) and we omit 
the proof which is much simpler. (4.14) is a bit different from (c). Integrate (7.4) over 
Vm+i to see that 

j y^+i\D,^^^pt{z,y)\i^{dy) 

< cr.,t-' liminf J y^^+,E,^^^ (^qsiym+u X^+'^m + z')'' + \No- z'|/^'])/^m+i(ciym+i) 
= cr.it-' liminf E,^^, ([(1 + z'y' + Wo - {{X^-^+'^y)) . 
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Now use the moment bounds in Lemma 3.3(d,e) to bound the above by 
ct-'^zr.,. ([(1 + ^T' + \No - z'\/z']{X^-^'^y). 



(7.34) 



The first term is trivially bounded by the required expression using Lemma 3.3 again. Using 
the joint density formula (Lemma 3.6), the Gamma power bounds (7.29), and arguing as 
in the proof of (b) above, the term in (7.34) involving A^o is at most 

dP-'^{z')-^¥. {\N-z'\{N\/ 1)P), (7.35) 

where N = N{z') is a Poisson random variable with mean z' . We have 

E {\N - z'\Nn^N>o)) < coiz' A (/)^/2+P) for aU ^' > and - l<p< 1/2. (7.36) 

For p <0 Lemma 3.3 of [BP] shows this (the uniformity for bounded p is again clear). For 
1/2 > p > use Cauchy-Schwarz to prove (7.36). Separating out the contribution from 
AT = 0, we see from (7.36) that (7.35) is at most 

ctP-\z')-^[e-'' z' + {z' A {z'y^''+P)] < ctP-\e-'' + 1 A {z'f-^/^) < ct^-^z' + Vf-^l'^. 



The result follows. 



□ 



Proof of Lemma 4.5(f). By the spatial homogeneity in the first m variables (7.26) we 
may use Lemma 4.4 to conclude 

\^Zi^lrr.^rVt{z,y)\ = I y D,.pt/2{x, -z^'^\ym+i)Dl^^^pt/2{-y^'^\ zm+i,x)i^{dx) 



Therefore 



/ 



\Dz^Dl^+^Pt{z,y)\i^{dz) 



< 




\D,^p,j2{x,-z^^\ym+iW'^^' 

X ^li+l\Dl^^^Pt/2{-y^'^\ Zm+ux)\n{dx)lXrn+l{dZm+l)- 

Use Lemma 4.2(a) to bound the first integral in square brackets and so bound the above 

by 

j zf,+,\Dl^^^Pt/2{-y^^\zm+l, X)\dx^^^ fXm+l{dZm+l) 

X {t + Xm+1 + ym+l)~^^'^qt/2{Xm+l, ym+l)lJ'm+l{dXm+l) ■ 
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The spatial homogeneity (7.26) imphes 

and so we conclude from the above that 
J ^m+i\Dz^Dl^+,Pt{z,y)\f^{dz) 

X {t + Xm+l + ym+l)~^^'^qt/2{Xm+l, Vm+l) IJ'm+l{dXm+l) 

m+l ; ym+ 

We have used Lemma 4.5(b) with g = in the next to last inequality and p < 3/2 in the 
last line. □ 

Proof of Lemma 4.5(e). For (4.18), use Lemma 4.4 and the spatial homogeneity (7.26) 
to bound the left-hand side of (4.18) by 

(m) 




X Zm+l\Dl^^^Pt/2{0,Zm+l,X^"''^ - Z^'^\xm+l)\l^idx) dz^'^^ ^rn+l{dym+l) ■ 

Use the substitution (for z'^^'>) w = x^'^'^ — z^"^^ and do the dx^'^^ firn+i{dym+i) integral 
first, using (4.13) to bound this integral by £4.5^^"^ (as p < 1/2). Now use (4.5) to bound 
the remaining dw fXrn+i{dxm+i) integral by C4,5t~^. 

The derivation of (4.19) is almost the same as above. One uses Lemma 4.2(b) now 
to bound the first integral. 

Proof of Lemma 4.5(d). The approach is similar to that in (b) as we integrate the 
bound in Lemma 7.1(a). There is some simplification now even with the higher derivative 
as g = 0. We use the notation from that proof, so that is the Gamma density in (7.27), 
Pniw) are the Poisson probabilities with mean to, and is a random variable with density 
Qn- Also let Bn be a Binomial (n, 1/2) random variable independent of F^. To ease the 
transition to Lemma 4.6 we replace t with s in this calculation. We also keep the notation 
z' = Zm+i/js, y' = yj^s. If 



01 (z', k) = (z' + 1)-Ml(fc<i)(l + (^0"') + kk>2)iz')-'[k +ik- z'/2)\ 
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and 



M^', k) = l(fc<2)(l + (zTn + kk>s)i^') [\k - + 3/c|/c - ^'1 + k], 
then by Lemma 7.1(a) and Fatou's lemma, the integral we need to bound is at most 

:= liminf cs-2[Ji(5) + J2{S)]. (7.37) 
By Lemma 3.6(b), 

Ji(5) < cs'/' / g5(2/m+i,2/) / E U 2-^n(^')</'i(^',fc)(^')'/'^^^V/^"'rfy 

poo ^ 

= cs3/2 / g,(y^+i,|/) Vp„(y')E(0i(r„i?„)r3/2)//7-irf^. (7.38) 

•^0 n=0 

The moment bounds (7.29) (the conditions there will be trivially satisfied now) give 



us 



E ((/>i(r,, 5jr3/2) < c(^p{Br, = o)E (r3/2) + p(b^ ^ i)e (r^^ + rl/^) 

+ E (i(B„>2)S„)E (r;^/2) + E (i(B„>2) (s, - r„/2)2r; 

+ i(„>2)E(E((s„ - r,/2)2|r,)r-3/2)). (7.39) 

The conditional expectation in the last term is [(F^ — n)^ + F„]/4. Therefore we may 
now use (7.31) and also (7.29) (as n > 2 and r = — 1/2 or —3/2 the conditions there are 
satisfied) to see that for n > 2 

E (E {{Br, - F„/2)2|F„)F-3/2) < (E ((F, - n)'T-^/^) + E (F" V2))/4 < en" V^. 
Insert this bound into (7.39) and conclude that 

E (01 (F,, Sn)r^/2) < c(2--(l + n5/2) + l(,>2)n-i/2) < ^ 
Therefore we can sum over n and integrate over y in (7.38) and obtain 

Ji(5)<cis^/^ (7.40) 
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where as always ci satisfies Convention 3.1. 

Lemma 3.6(a) and the argument leading to (7.38) shows that 



J2{S) < cs^/^ / 95(ym+i,y) J]Pn(y')E(02(r„,n)r3/2)yV7-irf^. (7.41) 

We have 

E((/>2(rn,n)r3/2) 

= 1(„<2)E ((1 + r;-)r3/2) + 1(^>3)E {r-^/\\n - r^l^ + 3n\n - r„| + n)). 

Some simple Gamma distribution calculations like those in the proof of (b), and which the 
reader can easily provide (recall Convention 3.1), show that the above is bounded by a 
constant depending only on Mq. As before by using this bound in (7.41) and integrating 
out n and y we arrive at 

J2{S) < c^s^'l''. (7.42) 
Insert the above bounds on Ji(5) into (7.37) to complete the proof. □ 



Proof of Lemma 4.6. Consider (4.22). The functions ^\ and ^2 are as in the previous 
argument. Argue just as in the derivation of (7.37) to bound the left-hand side of (4.22) 
by 



lim inf cs ^ 

5^0 



liminf cs-^[Ki(5) + K2{6)]. 



(7.43) 



Note we are integrating with respect to z^+i and not Hm+iidzm+i) as in the previous 
calculation. Lemma 3.6(b) implies that 



poo ^ 

Ki{6) <cs q6{ym+i,y)y2pn{y') 

•^0 n=0 



X {t/s) 



b/l 



t/s 



E{Mz',Bn))e- 



-dz'y^/^-^dy 



r(n + 6/7) 

<cs / q5{ym+i.y)Y.Pn{y') / E{c^i{z' ,Br,))-^-^dz'y'h-idy.(7.U) 
Jo -^0 [n + L) 



We have bounded r(n + 6/7) -"^ in a rather crude manner in the last line. 
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For < z' < t/s < 1 we have 

<c[p(5„ = 0)+P(S„ = l)(l + (z')"') 

+ 1(.>2)(^' + l)-\zT'^ (Sn + {Bn - z'/2f)] (/)"(n + 
< c[2--((/)" + n(l + {zT~') + lin>2){z'r-\n + {n - z'f){n + 1)-^ 

This, together with (7.44), shows that 

Ki{S) < ct. 
Next use Lemma 3.6(a) to see that 



< c. 



•^0 n=0 



X {t/s) 



t/s 





{z'Y 



r(n + 6/7) 



dz'y^/^-^dy 



<cs I q5{ym+i,y)^Pn{y') (f)2{z',n)-—-^dz'y''/''~^dy. 



n=0 



(n+1)^ 



As above, an elementary argument shows that for < z' < 1, (p2{z\n){z')^{n + 1) ^ is 
uniformly bounded in n, z' and also {b, 7) as in Convention 3.1. Hence, we may infer 

K2iS) < ct. 

Put the bounds on Ki{S) into (7.43) to complete the proof of (4.22). 

We omit the proof of (4.21) which is the first order analogue of (4.22) and is con- 
siderably easier. □ 

We need a probability estimate for Lemma 4.7. As usual X^'^'^^^ is the Feller 
branching diffusion with generator (3.1). 

Lemma 7.2. (a) F^iX^"^^^^ >w)< {wjzfl'^^ expj } for a\\w> z>^. 

(b) P^(Xi("'+^) <w)< {wjzfl'^^ exp| } for aU < «; < ^. 

Proof. This is a simple estimate using the Laplace transform in Lemma 3.3(c). Write Xi 
for Xf^^^^. If -{cit)-^ < A < 0, then 

P^(Xt >w)< e^"'E^(e-^^0 = e^^(l + At7)-''/^ expj \. 

l 1 + A^t J 
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If A > 0, then 

^z{Xt <w)< e^^'E.le-^^*) = e^^{l + At7)-^/^exp{^^^}. 

Now set A = ^^(^ in both cases. This is in (—('yt)~^, 0) if < 2; < and in [0, 00) if 
< z > w. A bit of algebra gives the bounds. □ 

Proof of Lemma 4.7. We will again integrate the bound in Lemma 7.2(b) over Zm+i 
but as g = we will use the function 

7/^(2;', n) = !(„<!) + l(„=i)(2;')-' + l(„>2)(2;')~^[(n - z'f + n]. 

We then have from Lemma 7.2(b), that for each z' > 0, 

= y J i^{ym+l<'W<Zm+l) "I" ^(.Zm+l<W<ym.+l))^rn+l 

X \Dl^+,Pt{z,y)\dz^'^^ Hm+iidym+i) 

< ct-^^^+iliminfE^^^, (^i;{z',No) J ^^(ym+i, ^i™^^^)(l(y^+i<«,<z^+i) 

+ '^izm+i<w<ym+i))f^m+l{dym+l) 

= ct-^2^+iE^^^, (^^/;(/,iVo)[l(«,<z^+i)l(^(-+i)<^) 

+ l(«,>a;^+i)l(x(-+i) >«,)]) 

< ct-2^j;+iE.^,,(V.(^',iVo)^)V2[l(„<,^,,)P,^,,(x(-+^) < wy/' 

In the third line we have used the a.s. and, hence weak, convergence of Xg^~^^^ to Xq^~^^^ 
as (5 — > and the fact that xj:"^~^^^ ^ w a.s. 
We have 



E {^P{z', Nof) < c((l + z')e-'' + iz')-^e-'' 

+ {z'r'i^ (l(Aro>2) (A^o - z')') + E (l(^„>2)Ar2)]) . 



An elementary calculation (consider small z' and large z' separately) shows that the term 
in square brackets is at most c{z')'^. Therefore we deduce that 

E{i;{z',Nofy/^<ciz')-'- (7.46) 
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If we set w' = w/^t, then this, together with (7.45) and Lemma 7.2 allows us to conclude 
that 

J J i^m+l) fJ-m+lidZm+l) 

<ct J z^^^-^ ' \z) {w/zm+i) ' ^exp(^ ^ — jdZm+1 

= ct^-^+'l\w'f/^^K^_2+^,/^,{w'). (7.47) 

A simple calculation using the obvious substitution x = {y/z' — s/w')'^ shows that for any 
£ > there is a co{e) such that 

Kr{w') < co[l(^/<i) + (u;')'^+^/^l(«,'>i)] for aU > and - 1 + £ < r < 

Our bounds on p and Convention 3.1 imply that r = p-2 + 3b/ 4^ e [-1 + 3(4M^)-^ M^] . 
Therefore the left-hand side of (7.47) is at most 

□ 

8. A Remark on the Two-dimensional Case. 

As has already been noted, the proof of Proposition 2.2 (by far the most challenging 
step) simplifies substantially if d = 2. As this is the case required in [DGHSS], we now 
describe this simplification in a bit more detail. 

Recall the three cases (i)-(iii) for d = 2 listed following Theorem 1.4. As noted there, 
the case £^ = is covered by Theorem A of [BP] (with d= 2) without removing (0, 0) from 
the state space, so we will focus mainly on the other two cases here (but see the last 
paragraph below). In these cases the localization in Theorem 2.1 reduces the problem to 
the study of the martingale problem for a perturbation of the constant coefficient operator 

2 

A' = Y.b^D,,+^^x,Dl, (8.1) 

i=l 

with resolvent Rx and semigroup Pt- Our job is to establish Proposition 2.2 for this 
resolvent. 
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For f eVq, we have 



M°i?A/||2 = ||Ai?A/-/||2<2||/||2, 



(8.2) 



the latter by CoroUary 2.12. We may therefore remove the term X2{R\f)22 from the 
summation in (2.6) because L^-boundedness of this term will follow from the other three 
and (8.2). Recall the required boundedness of any of the terms was reduced to (2.28) by 
Cotlar's Lemma. (Proposition 2.10 was only used to ensure the operator Tt was bounded 
on so that Cotlar's Lemma may be employed in the proof of Proposition 2.2. This 
boundedness, as well as the bound of ct~^, is also implied by (2.28) with s = t and 
the elementary Lemma 3.16. So we only need consider (2.28).) For the two derivatives 
involving xi, (2.28) was fairly easily checked in Lemma 4.3 thanks to the "explicit" formulae 
(4.7) and (4.8), and the bound in Lemma 3.3(f). 

It remains only to check (2.28) for ■ This was done in Proposition 4.9, using 
only Lemma 4.6(a) and in fact only used (4.14) for p = and (4.13) for p <0. These proofs 
in turn were fairly simple consequences of part (c) of the key Lemma 7.1. (Admittedly the 
proof of (4.13) was omitted, being much simpler than that of (4.15).) As (7.4) was a trivial 
consequence of (7.3) (recall (7.26)), we have essentially reduced the two-dimensional case 
to the proof of (7.3). To justify our earlier statements, that this really is much simpler 
than that of (7.1), we give the proof. At the risk of slightly lengthening the argument we 
will take this opportunity to explicitly write an integration by parts formula which was 
implicit (and hidden) in the more complicated setting of Lemma 7.1. Recall that m — 1 
(the proofs below are the same for general m). It — xi^"* ds, and (see (3.26)) 



A^o = Ao(t) is the Poisson variable in Lemma 3.4. 

Proposition 8.1 (Integration by Parts Formula). If / : R+ x R — > R is bounded and 
Borel, then 




{No - {Z2M) 



Gt,,J{IuXi^^)')+E^{t,z,f), (8.3) 



where Ei is given by (8.6) below and satisfies 



E^{t,z,f)\<E 



(/ / y"|/(a;i,xf))|4/rV470(,+,^)C 



xi — zi — hit) 




(8.4) 
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Proof. By (3.25) 

+ E (y [Gt,,, f ( J XP + usds, ) + ut) 



= E^{t,z,f) + E2{t,z,f). 



(8.5) 



Now use 



dpt 
dt 



and (by some calculus) 



[z) = {zH-'-l){2t)-'pt{z), 



dpt 4 



together with the Fundamental Theorem of Calculus, to obtain 

■(xi-^i-6?t)2 



E^{t,zJ) =E 




(2(« + ^t)) 



-1 



X 



27?(« + /t) 

P270(u+/t)(a;i -^1 -^>iOl(„<j*^^d^)C^a;icZwl(^^=o)t^No(z^)), (8.6) 



and 



\E^{t,z,f)\<E 



Z2 




/(a;i,xf))|4(« + /t)-V47?(u+/*)(a^i --^1 -^?*) 



The latter inequality gives (8.4). 

For E2 we use the decomposition in Lemma 3.4 with p = 0. Sn and i?„ are the sum 
of the first n of the and r^, respectively, and we continue to write Pn{w) = e~'^w'^/n\ 
and z' = zij^t. Then Lemma 3.4 allows us to write 

E2(t,^,/) = No(i^t > 0)E 

- Gt,,J{h{t) + RNo,X'o{t) + Smo)) 

00 

= (7^)"' Y.(Pn-l{z') -pn{z'))E,^{Gt,,J{h{t) + Rn.X'^it) + S^)) 



n=0 

00 



= (7*)-' - ^')(^')"'E,,(G,,,,/(J2(t) + Rn,X'^{t) + 5J) 

n=0 

= (7t)-^E,,((iVo - ^')(^')-'G',,,,/(/„xf ))). (8.7) 
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In the last line we have again used Lemma 3.4 to reconstruct X^'^\ (8.7) and (8.5) complete 
the proof. □ 



Remark 8.2. Since \Gt,zif \ < 
\D.,Ptf{z)\ 

\Nq - z' 



the above implies the sup norm bound 



<(7t)-'E,,( 



z' 



We have used Lemma 3.3(f) and (3.17) in the last. This gives a derivation of (4.4). More 
importantly we can use the above to derive an bound which will allow us to take / = 5y. 
Recall that qt{x,y) is the transition density of X*^^^ with respect to y^/'^~^dy. 



Corollary 8.3. If / : R+ x ] 

\D,,Ptf{z)\dzi 



is bounded and Borel, then 



1/(^1, xf))|cZ^i 



\NQ-Z2ht\ ^ |^£2 ^^^-ijy 



It 



Proof. Note first that / X)|(izi < / \ f{xi,X)\dxi^ and then integrate over zi 

in Proposition 7.1 to see that the above integral is at most {z' = Z2/^t as usual) 

|iVo-^'| 



Z2 



+ E 



22 



|/(xi,xf))|c/a; 



\f{x,,xi^^)\dxi4Ii' 



ds d'H\ 



o(^)). 



(8.8) 



Use Lemma 3.3(f) and (3.17) again to bound the last term by 

4c3.2(t + Z2)-^¥.,^ ( j \f{xu ))|dXi ) . 

Use this in (8.8) to derive the required bound. 



□ 



Proof of (7.3). Let p'^{zi,Z2) = P5{zi —2/1)175(2/2,-22) (bounded in z by Lemma 3.3(a)). 
Then (3.30) shows that lim^^o ^z2^t/^'^(^) = D^^pt {z,y). Apply Fatou's Lemma and 
Corollary 8.3 to conclude 



/ 



\Dz2Pt{z,y)\dzi 



< liminf C7 3t"^E 



\fy'\z„xn\dzi 



(2)^ 



|iVo - ^2/7^1 



22 /t^ ' ^7* 
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The required result follows. 



□ 



If we wanted to include the case £^ = to make the above "short proof self- 
contained, then we need to consider Proposition 2.2 and hence (4.1) and (4.2) for the 
case 

1=1 

The associated semigroup Pt = 11^=1 Qt is a product of one-dimensional Feller branching 
(with immigration) semigroups with transition densities given by (3.3). As in the the last 
part of the proof of Proposition 2.14 at the end of Section 4, (4.1) and (4.2) reduce easily 
to checking (4.1) and (4.2) for each one dimensional Q]. In the first part of the proof 
of Proposition 2.14 (in Section 4) we saw that these easily followed for each differential 
operator by projecting down the corresponding result for (as in (8.1)) to the second 
coordinate. This was checked in the "short" proof above for the the first order operators. It 
therefore only remains to check (4.1) and (4.2) for Dx = xD"^ and qt in place oipt- As in the 
proof of Proposition 4.12, we must verify (4.33), (4.34), (4.24), and (4.25) for this operator 
and one-dimensional density. These, however, can be done by direct calculation using the 
series expansion (3.3)-the arguments are much simpler and involve direct summation by 
parts with Poisson probabilities and elementary Poisson bounds. 
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